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ABSTRACT 


The  determinantal  equation  for  a  narrow  tape  helix  is  derived 
of  two  methods,  and  complex-valued  solutions  for  the  phase  constant 
are  obtained. 

The  complete  k  -  B  diagram  (Brillouin  diagram)  is  given  as  a 
function  of  tape  width  and  pitch  angle.  In  order  that  the  solutions 
be  continuous  functions  of  k  and  B,  it  is  necessary  to  change  branches 
of  the  square  root  which  appears  in  the  determinantal  equation.  A 
discussion  of  solutions  which  are  physically  admissible  as  complex 
wave  solutions  is  given,  and  the  phase  constants  corresponding  to  the 
complex  wave  solutions  are  used  to  represent  the  current  on  a  helix. 

Two  source  problems  are  investigated,  one  an  infinite  helix,  the  other 
a  finite  helix.  Comparison  with  experiments  is  made  with  good  agreement. 
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1.  INTRODUCTION 


Analysis  for  the  helix  was  first  reported  in  1897  by  Pocklington''"  who  as¬ 
sumed  a  very  thin  perfectly  conducting  wire  and  formulated  an  integral  equation 
whose  approximate  solutions  predicted  an  axial  phase  velocity  equal  to  the  veloc¬ 
ity  of  light  c  for  very  low  frequencies  and  c  sin  t ,  where  41  is  the  pitch  angle 
of  the  helix,  for  larger  frequencies.  These  results  agree  with  the  results  re¬ 
ported  in  the  present  work  for  the  limiting  cases  of  very  narrow  tape  and  low  fre¬ 
quencies  . 

2 

In  1910  Nicholson  formulated  the  helix  problem  exactly  but  was  forced  to 
make  unreasonable  approximations  to  obtain  a  solution. 

3 

Ollendorf  in  1926  reported  the  solution  of  the  sheath  helix  for  the  mode  in 

which  no  angular  variation  of  the  fields  exist.  Solutions  for  the  higher  order 

4 

modes  of  the  sheath  helix  were  given  by  Phillips  and  Malm  . 

A  more  complete  analysis  which  was  motivated  by  the  use  of  the  helix  in 

5 

traveling  wave  tubes  and  as  antennas  was  given  by  Kornhauser  in  1949. 

An  excellent  history  and  an  extensive  bibliography  through  1955  is  given  by 

6  7 

Sensiper  who  also  explains  the  relationship  of  his  thesis  to  the  total  litera¬ 
ture.  In  these  publications  he  formulated  and  solved  the  determinant a 1  equation 
for  the  tape  helix  for  the  phase  constant  13.  He  limited  himself  to  the  real¬ 
valued  solutions  of  this  equation. 
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Following  Sensiper,  Pierce  and  Tien  obtained  an  approximate  solution  of  the 

9 

determinantal  equation  of  the  tape  helix  using  Pierce's  coupled  mode  theory. 

Abstracts  of  66  papers  and  reports  on  the  helical  beam  antenna  through  May 
1959  are  given  by  Wong  and  Thomas'^. 

None  of  the  works  referred  to  above  indicated  the  possibility  of  the  exis¬ 
tence  of  complex-valued  solutions  of  the  helix  determinantal  equation.  In  fact, 
Sensiper  and  also  Pierce  voiced  doubts  as  to  the  existence  of  complex-valued 


2 


phase  constants,  Senslper' s  proof  on  the  non-existence  of  complex  solutions  was 
based  on  the  condition  of  total  finite  energy. 

The  existence  of  the  complex  roots  and  their  importance  to  the  helical  antenna 

problem  was  first  pointed  out  by  Mittra^\  In  that  paper  he  also  explained  as 

a  consequence  of  the  presence  of  a  particular  symmetry  m  the  helix  structure  which 

makes  its  effective  period  to  be  infinitesimally  small,  why  the  determinantal 

equation  of  the  helix,  unlike  the  equation  for  other  periodic  structures,  does 

not  have  periodic  solutions.  As  a  further  consequence  of  the  above  geometry  it 

was  shown  that  there  exist  no  harmonic  terms  in  the  representation  of  the  current 

i  P  ws 

along  the  tape  and  that  a  single  term  of  the  type  e  (s  -  distance  along  the 

tape,  Pw  =  phase  constant)  is  sufficient  to  represent  this  current. 

A  major  contribution  of  the  present  work  is  the  detailed  study  of  the  deter¬ 
mination  and  application  of  the  complex  solutions  of  the  helix  determinantal 
equation.  Except  for  a  brief  mention  by  Mi+tra,  a  discussion  of  these  solutions 
has  not  been  reported  elsewhere.  The  complex  valued  solutions  of  the  helix 
equation  have  been  found  very  useful  in  explaining  the  current  distribution 
on  the  helix,  both  uniform  and  log-periodic,  and  for  predicting  the  radiation 
pattern  of  a  helical  antenna.  A  further  discussion  on  the  interpretation  and 
usefulness  of  the  complex  solutions  of  the  determinantal  equation  in  an  open 
periodic  structure  appears  later. 

An  alternative  approach  to  the  solution  for  the  current  distribution  on  the 
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helical  structure  has  been  reported  by  Patton  .  He  obtained  the  Fourier  trans¬ 
form  of  the  current  distribution  on  the  semi-infinite  bifilar  helix  by  a  Wiener- 
Hopf  technique.  The  Fourier  transform  of  the  current  was  related  to  the  radi¬ 
ation  pattern,  and  the  calculated  patterns  were  compared  with  experimental,  results 
also  obtained  by  Patton,  The  effect  of  wire  size  was  studied. 

13  14  15 

Three  related  papers  have  recently  been  written  by  Maclean  ’  ’  .  In  the 


first  he  compares  the  different  approaches  to  the  helical  antenna.  The  results 


3 


of  the  report  state  that  a  "simple  engineering  approach"  ('i.e.,  one  without 

Bessel  Functions)  "is  adequate  for  mosi  purposes/'  The  second  paper  is  a 

theoretical  study  assuming  that  waves  travel  both  axially  and  helically  along 

the  antenna.  It  is  shown  that  the  existence  of  the  helical  waves  is  responsible 

for  what  are  the  cut  off  frequencies  of  the  antenna.  Study  is  also  made  of 

the  effect  of  the  ground  plane.  The  last  is  an  analysis  based  on  the  sheath- 

helix  mode.  He  obtains,  among  other  things,  the  determinantal  equation  for 

the  sheath  helix  with  a  concentric  perfectly  conducting  core.  The  results  are 

given  for  relative  phase  velocity  when  the  fields  have  one  axial  variation, 
~in<j> 

i.e.,e  ,  n  =  1 . 
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Kraus  in  1948  and  later  Kornhauser  calculated  the  radiation  pattern  of 
a  helical  antenna  operation  on  the  axial  or  beam  mode  by  assuming  a  sinusoidal 
current  with  an  empirically  determined  phase  constant  which  corresponded  to  a 
slow  wave  on  the  antenna.  They  obtained  satisfactory  results  by  using  the  em¬ 
pirically  determined  phase  constant  which  may  now  be  obtained  by  using  rhe  phase 
constant  determined  from  the  solution  to  the  de1  e rminantal  equation. 

18 

The  current  distribution  on  a  helical  antenna  was  measured  by  Marsh,'  and 

he  was  able  to  fit  his  measured  results  well  by  assuming  there  existed,  two  waves 

on  the  antenna,  one  corresponding  to  a  complex- valued  phase  constant  with  the 

real  part  representing  a  wave  traveling  at  the  velocity  of  light  and  the  other 

corresponding  to  a  real  valued  phase  constant  representing  a  slow  wave.  His 

empirically  determined  phase  constants  agree  well  with  the  complex- valued  phase 

constants  which  are  solutions  to  the  determinanta.l  equation. 

Now  return  to  a  brief  discussion  of  complex  solutions.  The  solutions 

corresponding  to  complex-valued  phase  constants  are  not  proper  modes  since  the 

fields  corresponding  to  these  modes  do  not  satisfy  the  radiation  condition 

19 

everywhere  a+  infinity.  However,  Warcuvitz  suggested  that  these  solutions 


might  be  used,  to  approximately  represent:  the  near  fields,  and  used,  as  an  alter- 
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nate  to  the  continuous  eigen-value  spectrum  and  as  an  approximate  representation 
of  the  far  fields  of.  a  structure  in  a  certain  range  of  observation  points.  The 
solutions  corresponding  to  the  improper  modes  were  given  the  name  leaky  modes  by 
Marcuvitz  who  showed  that  the  zero  of  the  de'etnunant ai  equation  which  corres¬ 
ponds  to  the  leaky  mode  is  on  the  improper  Riemann  sheet  in  t  tie  phase  constant, 
plane.  The  far  fields  may  also  be  calculated  by  Kirchboff  type  Integra' ion 
by  treating  the  near  field  representation  in  terms  of  the  leaky  wave  fields 
as  equivalent  sources. 

Radiation  properties  of  a  class  of  periodic  structures  may  be  conveniently 

analyzed  in  terms  of  their  Brillouin  (k-  8)  diagram  provided  this  diagram 

exhibits  both  real  and  complex  solutions  in  slow  as  we  1 1  as  fast  wave  regions. 

Early  workers  on  periodic  structures  were  primarily  interested  in  'he  surface 

wave,  (real  valued  solution  to  the  determi nani al  equation),  on  open  periodic 

structures,  and  it  was  customary  to  name  the  fast  wave  regions  forbidden,  regions. 

Recently  it  has  been  well  demonstrated  that  not  only  do  there  exist  complex-valued 

solutions  for  the  phase  constant  of  some  structure  in  the  forbidden  region  where 

8  <  k  but  also  in  the  slow  wave  regions  where  6  >  k.  Since  leaky  waves  have 

r  r 

become  synonymous  with  fast,  waves,  it  is  suggested  here  that  the  name  complex 
waves  be  associated  with  the  waves  corresponding  to  complex-valued  phase  con¬ 
stants  in  fast  or  slow  wave  regions  o.r  that  the  term  leaky  wave  be.  generalized 

to  include  both  slow  and  fast  waves.  The  physical  interpretation  of  complex 

19  20 

waves  was  given  by  Marcuvitz  .  Olner  gives  an  excellent.  discussion  on  the 
range  of  permitted  observation  points  and  me  position  of  the  complex  wave 
zero  of  the  determinantal  equation.  In  a,  class  of  structures  excited  by  a 
given  source,  where  the  dominant  pact,  of  the  field  has  its  total  representation 
in  terms  of  the  continuous  spectrum  of  spatial  frequencies,  complex  waves 
have  been  found  useful  in  representing  the  fields.  This  alternate  repre- 
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sentation  is  approximate. 

2 1  22 

Goldstone  and  Oliner  •’  obtained  good  results  using  complex  waves 

in  the  analysis  of  a  class  of  leaky-wave  structures.  'The  relationship 

between  complex  waves  and  the  radiation  pattern  of  leaky  wave  antennas  is 
23 

discussed  by  Hessel.  Oliner  also  gives  other  examples  of  the  use  of  com¬ 
plex  waves  in  electromagnetic  phenomena. 

To  date,  very  few  structures  have  been  analyzed  for  their  complete 
Brillouin,  (k-B),  diagram,  other  than  on  an  approximate  and  perturbat i onal 
basis.  This  generalization  is  excepting  some  idealized  structures  such  as 
the  constant  reactance  or  modulated  reactance  surfaces,  the  anisotropic  sheet 
structures  and  a  few  others.  Although  the  analysis  of  the  above  structures 
is  quite  useful,  their  application  to  practical  devices  is  good  only  in  an 
approximate  sense. 

The  tape  helix,  which  is  very  much  a  practical  structure,  yielded  to 
solution  for  its  complete  k-13  diagram  which,  with  its  applications  to  heli¬ 
cal  antennas,  is  the  subject  of  this  report.  The  determinantal  equation  for 
a  narrow  tape  helix  may  be  formulated  in  at  least  two  different  ways.  One 
involves  the  Floquet  expansion  of  the  fields  and.  the  assumed  current,  with 
unknown  phase  constant,  on  matching  boundary  conditions  at  the  cylindrical 
surface  which  contains  the  tape.  The  second  method  is  based  on  an  integral 

equation  formulated  for  the  current,  distribution  on  the  tape.  The  second 
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method  is  similar  to  the  method  used  by  Kogan  on  the  round  wire  version  of 
the  helix.  The  same  determinantal  equation  is  obtained  by  both  methods.  The 
determinantal  equation  and  its  solutions  for  the  phase  constant  are  ex¬ 
tensively  studied  for  real-valued  as  well  as  complex-valued  phase  constants. 

The  solutions  in  the  k-6  plane  are  continuous.  In  order  that  the 
solutions  be  continuous  it  is  necessary  to  change  branches  of  the  square 
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root  in  the  argument  of  the  IK  product  in  determinantal  equarion.  Different 

branches  give  .rise  to  different  solutions  which  are  physically  interpreted. 

The  root-tracking  procedure  on  different  Riemann  sheets  which  is  outlined 

here  may  find  applications  in  the  study  of  other  structures  where  the  locus 

in  the  k-S  plane  appears  to  be  discontinuous. 

The  complex-valued  solutions  to  the  determinantal  equation  are  used  in 

the  analysis  of  two  different  source  problems.  The  first  source  problem 

obtains  the  solution  for  the  relative  amplitude  of  the  waves  corresponding 

to  real  and  complex  roots  for  an  infinite  helix  excited  a+  the  center.  The 

second  considers  the  same  problem  for  a  finite  helix  excited  at  the  center. 

A  variational  principle  is  used  to  formulate  the  latter  source  problem. 

Theoretically  computed  values  of  the  real  and  complex  solutions  are  used  to 

18 

explain  the  experimentally  measured  current  distributions  by  Marsh  .  As 
previously  mentioned,  complex  wave  analysis  given  in  this  report  may  be  used 
to  interpret  Marsh's  results  which  were  empirical. 

It  is  felt  that  the  current  distribution  on  the  tapered  version  of  the 
helix,  viz.,  the  equiangular  spiral  on  a  cone,  may  also  be  explained  by  ex¬ 
tending  the  present  work  to  tapered  structures. 
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2.  THE  DETERMI NANTAL  EQUATION 


2.1  Introduction 

In  this  chapter  the  determinantal  equation  for  the  infinite  tape 
helix  will  be  derived  by  two  methods.  The  functions  of  a  complex 
variable  which  appear  will  be  multivalued,  and  determination  of  the 
appropriate  branches  will  be  discussed.  The  geometry  of  the  tape  helix 
used  is  shown  in  Figure  2.1. 

The  tape  lies  in  a  cylindrical  surface  whose  radius  is  "a"  and  whose 
axis  is  the  z  axis.  Perfect  conductivity  and  infinitesimal  radial 
thinness  are  assumed;  however,  the  width  of  the  tape  is  taken  as  d  as 
measured  in  the  z  direction.  Points  on  the  surface  of  the  cylinder  are 
given  by  the  coordinates  and  The  azimuthal  angle,  4>,  locates  a 
point  on  the  center  line  of  the  tape  as  given  by  the  parametric 
equations 


x  =  a  cos  4> 
y  =  a  sin  <}> 


P  <f> 


where  p  is  the  pitch  of  the  helix.  The  coordinate  £  is  the  distance 
in  the  z  direction  from  the  point  determined  on  the  center  line  of  the 
tape  by  specifying  <j>. 

The  equation  of  the  tape  may  then  be  written  as 


-oo  <  4“  <  °° 
6  6 
-2<  ^<2 
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2.2  Derivation  by  Application  of  Boundary  Conditions 

The  determinantal  equation  for  the  tape  helix  is  derived  in  this 

25 

section  by  expanding  the  fields  and  current  in  space  harmonics 

and  then  matching  boundary  conditions  at  r  ■=  a. 

The  helix  shows  a  great  amount  of  symmetry.  If  one  displaces  the 

infinite  helix  a  distance  p  in  the  z  direction,  the  helix  transforms 

into  itself.  The  translational  period  is  p.  If  the  helix  is  rotated 

through  an  angle,  <j>  =  2 TT}  the  helix  transforms  into  itself.  The 

rotational  period  is  2ir.  These  symmetries  are  not  the  only  symmetries. 

In  addition,  any  arbitrary  translation  z  followed  by  a  rotation  <j*  where 
z 

4>  =  —  transforms  the  helix  into  itself.  Therefore  the  period  in 

o 

p  _  1  1  z  -  1 

4  =  z  -  p4>  ,  where  =  —  ,  is  zero..  Now  if  C  =  z  -  p  <f>  and  is  fixed 

then  ^  is  periodic  with  period  p.  If  z  is  fixed  ^  is  also  periodic  with 

period  p.  This  property,  the  generalized  period  being  zero,  is  a  unique 

property  of  the  helix.  Of  course,  uniform  structures  have  zero  translational 

period.  This  unique  symmetry  property  of  the  helix  has  important 

consequences  which  will  be  discussed  later. 

The  field  equations  used  are  expressed  in  circularly  cylindrical 

coordinates.  Starting  from  the  source  free  Maxwell's  equation  one 

obtains  the  wave  equation  for  E  and  H  as 

z  z 


I  9 

r  <5"r 


34'.  1 

(r  3--)  t  ~ 
r 


where  4*  =  E  or  H  and  is  a  function  of  r ,  4>,  z  and  t. 

7.  7.  J  ■ 


Assuming  the  time  variation  as  sinusoidal  and  that  Equation  (2.1)  may 
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be  separated  as 


(r,  <j>.  z,  t)  =  R(r)  |>(<i>)  Z(z)  e1Ut 


one  obtains 


where 


d2|  2  - 
— -  +  n  4>  =  0 

d4> 


and 


k2  =  u2jUf. 


A  solution  to  Equation(2„3) is 


A1e' 


i  n4> 


and,  since  the  fields  are  periodic  in  <j>  with  period  2T r,  n  must  be  an 
A  solution  of  Equation(2 . 4) is 


Z  =• 


Bm 


.11 

A  e 


z 


(2,2) 


(2,3) 


(2,4) 


(2,5a) 


integer. 


(2 „ 5b) 
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The  solution  for  Equation(2.2)  is 


where 


(2,5c) 


T 

m 


2 

(k  p) 


If  the  fields  are  to  be  finite  at  r  =  0  or  r  =  oo,  then  the  fields 
are  of  the  form 


^(r^  z) 


oo 

=  £ 

m,  n=-oo 


(2.6) 
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Floquet's  Theorem  states  that  the  fields  at  z  =  z  +  p  are  related 

to  the  fields  at  z  by  a  complex  constant.  If  B  =  B  4.  m  then  Equation  (2 .6) 

m 

may  be  written 


-iB 


oo 


4,(r,4';Z)  =  e  P  E 


m ,  n=-oo 


mn 


In(T*  5) 

K"  C' ;) 


-in4>  p 

e  e 


(2.6a) 


Now  Equation  (2.6a)  satisfies  Floquet’s  Theorem  since 


^(r,  <j>,  z  +  p)  =  e  +(r,<j^z) 


When  z  is  changed  to  z  +  z  and  <)>  changed  simultaneously  to  6  +  cfa  where 

o  o 
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z  =  p  4>  ,  the  fields  must  be  the  same  except  for  a  complex  constant, 

o  o 

The  property  that  the  fields  differ  only  bv  a  complex  constant  after 

the  above  combination  of  translation  and  rotation  is  the  result'  of 

26 

"The  Generalized  Floquet  Theorem  In  "The  Generalized  Floquet 

Theorem"  use  is  made  of  the  periodicity  not  only  in  the  translation  alone 
or  in  the  rotation  alone  but  in  the  combination  of  translation,  rotation, 
and  other  kinds  of  reflection  symmetries. 

The  Generalized  Floquet  Theorem  is  satisfied  for  the  helix  whose  fields 
are  described  by  Equation  (2.6a)  if  and  only  if  n  =  -m.  Equation  (2.6a) 
becomes 

S 

-1  —  z 

oo 

+  (r,<j>,z)  =  e  P  £  A 

n 

n-~oo 


Let  the  logitudinal  fields  be 


where  i  and  0  refer  to  r  <  a  and  r  >  a  respectively. 

By  applying  Maxwell's  equation  the  other  fields  are 


(2.8) 


(2.9) 


4* 
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The  boundary  conditions  for  the  tape  helix  at  r  =  a  are-, 

1)  Tangential  electric  field  is  continuous  everywhere  on  the  cylinder 
r  =  a . 


I 
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2)  Tangential  magnetic  intensity  is  discontinuous  by  the  surface 
current . 

The  surface  current  will  be  assumed  to  hate  the  following  form 


J 


11 


J  -64 

Oc* 


-  <j>  <-  z  <•-<)>  j.  5 

p  p 


J 


1 


=  o 


elsewhere  in  0  <  z  <,  p 


(2.14) 


Such  an  assumption  is  justified  because  the  fields  at  4  -s-  4,  are  related 

to  the  fields  at  4  exactly  in  the  manner  in  which  the  fields  at  4*  4  4-^ 

are  related  to  the  fields  at  cj)1 ,  where  4,  4*  and  4j  are  arbitrary.  The 

reason  for  the  previous  statement  is  that  if  the  helix  is  infinite  then 

one  cannot  distinguish  between  the  fields  at  4  and  4’*'  except  by  the  complex 

constant  as  shown  in  Equat ion (2 . 7 0  Later.,  in  Chapter  5,  it  will  be  shown 

that  the  solution  for  an  infinite  helix  fed  at  the  origin  will  be  of  the 
-1  31  4 

form  e  .  This  current  is  confined  to  the  tape  and  is  flowing 

only  in  the  parallel  direction.  In  addition  the  current  has  a  phase 
shift  only  when  4  of  the  center  line  is  changed.  Thus  the  current 
satisfies  Floquet's  Theorem.  If  one  is  to  be  successful  in  matching 
boundary  conditions  the  current  must  be  expanded  in  the  form  of  Equation  (2.7), 
viz. 


oo 


Ju  =  e  P  L  *  jlln  e  p  (2.15) 

n=-oo 


This  form  for  the  current  is  identical  to  the  form  for  the  fields, 


15 


Equation  (2.7),  and  permits  the  matching  of  the  boundary  conditions. 

If  one  examines  Equat  ion  (2 . 15),  Equation  (2 . 14)  would  be  rewritten  as 


R  z  .  B  r 

-i  Is  -  ~ 

Je  p  e  p  0  <  ^  <  S 

o 

0  elsewhere 


(2.16) 


The  current  of  Equation  (2 . 16) l s  easily  expressed  in  a  Fourier  series 
of  the  form  of  Equation  (2  15).  The  coefficients  are  found  to  be 


jlln  =  -  f  J 
p  J  o 


B  *  n  r 
-1  —  b  1  -  b 


e  p  e  p  dr 


6  6  ^ 

s  -l  B  -  sin  - 

J  0  n  - 

,  o  2p  2p 

Jlln  =  —  e  “  g  5— 

n 

2p 


(2.17a) 


j4*n  =  jlln  cos  + 


(2.17b) 


=  j 


=  jlln  sin  41 


(2.17c) 


Now  since  all  the  fields  and  currents  are  expanded  in  similar  series  of 
linear  independent  functions,  the  boundary  conditions  must  apply  to  each 
term  of  the  series  or  space  harmonic.  The  boundary  conditions  are 
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,  /  .  nSn/-  .  .  .  / 

A  °K  1  (t  -A  -  - - E—  B  °K  /r  -  A  *1  1  [T  - 

n  n  \  n  p/  _n  n  n  \  n  P/  n  n  \  n  p/ 


+  n  “  B  1I  (t  -  jlln  sin 

__p - *  n  l n  p/ 


(2  19d) 


The  above  four  equations  can  be  solved  for  A  1,0  and  B  1,C>.  These  solutions 

n  n 

yield 


K  T 
n  n 


.  T  a 
a\  n 


P/  P 
iu* 


n  di 

—  sin  ^ 


(2.20a.) 


3  1  =  -  K  1  (r  i\  I*-* 
n  n  I  n  - 

\  P/  P 


jlln  cos  ^  (2.20b) 


I  T  - 
n  n  - 

A  °  =  A  1  _ V-  .  P 


(2.20c) 


If  a' 
I  r  - 

U  °  a  1  n  n  - 
B  =  B  \  p, 

n  n  — ; — '• - 


(■•  I 


(2 . 20d) 


where  use  of  the  Wronskian.  I  '''(x)  K  (x)  -  I  (x)  K  ''"(x)  ■=  was  made. 

;  n  n  n  n  ^ 
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The  component  of  electric  field  parallel  to  the  helix  may  now  be 
written  as 


11 


oo 

=  L 
n=-» 


1  In 


(2 „ 21) 


where 


E  =  E  sin  4*  +  E  ,  cos  41  (2.22) 

lln  zn  q>n 


If  Equations  (2 . 8),  (2  .  id)  and  (2 . 20)  are  substituted  into  Equation  (2 . 22)^^^ 
is  obtained  viz. 


If  Equation  (2 . 23)  is  substituted  in  Equation  (2 . 2l)  and  E  set  equal  to 
zero  along  the  center  line  of  the  tape,  i.e.,  !»  =  0,  then  the  determinantal 
equation  for  6  is  obtained 
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7 

Equation  (2.24)  has  been  found  by  Sensiper  .  It  was  reproduced  here 

using  the  concept  of  generalized  periodicity.  The  factor  outside  the 
summation  cannot  be  zero.  Therefore  both  sides  of  Equation C2 .24)  may  be 
divided  by  the  non-zero  factor.  If  the  following  relations  among  the 
Bessel  functions 


I  (x)  K  (x) 
n  n 


■I  K  - 
n-1  nil 


1  .-,K  1 

ni-1  n-1 


I  1(x)  K  1(x) 
n  n 


I 

n+ 


1 


K 

n+ 


1 


are  substituted  into  Equation  ^ .  2^1,  one  obtains 


0  =  L 


~~  -  k2)  a2  I  K  -  -  ctn2  +  (i  K  ~I.  K  \  D 

-2  /  n  n  2  ^  n-1  n-1  n+1  ni-lj  n 


(2.25) 


where 


D 

n 


sin 


B  6 

n 


_ 2  p 

B 

n 


2p 


The  argument  of  each  IK  product  in  Equations  (23 ,  24  and  25)  is 


If 


A  . 

l 


oo 

L  I  (T  -)  K  (T  --)  D 
n+i  v  n  -  '  n+i  n  -  n 
n=~oo  p  p 


then  Equation (2. 25)  may  be  written  as 


2  -2 
k  P 


sin 


k2p2ctn2  ^ 


A  +  A  -  2A 
+  1  -1 _ o 

2A 


(2.26) 


Equation(2 . 26) is  the  determinantal  equation  for  the  tape  helix  derived  by 
assuming  a  surface  current  uniform  across  the  tape  with  phase  shift  along 
the  length  of  the  tape.  The  form  of  Equation (2 .26) is  suitable  for  studying 
the  behavior  of  B  as  a  function  of  k. 

2,3  Derivation  by  Transform  Techniques 

In  this  section  the  integral  equation  for  the  electric  field  on  an 
infinite  helix  is  derived.  The  determinantal  equation  is  then  obtained 
by  taking  the  Fourier  transform  of  the  electric  field  along  the  helix. 
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A  current  I(4>)  is  assumed  along  the  center  line  of  the  helix  and  then 
the  electric  field  is  calculated  by  finding  first  the  vector  potential 


A 


(2.27) 


and  second  the  electric  field 


E  = 


tuA  + 


grad  div  A 
1  uP* 


(2.28) 


The  calculation  gives 


E 


I 


-ikR 

e 

R  ~ 


dL 


+ 


-ikR  I- 

/v  A  6  . 

a  a  — - —  1 

R  R  R  L 


-f-  + 

K  R 


dL 


(2.29) 


where  'a  is  unit  vector  pointing  toward  the  observation  point  from  the 

R 

source  point.  The  component  of  electric  field  parallel  to  the  center 
line  is  given  by  E  -  E  '3^.  On  the  surface  of  the  cylinder 


-Ujfi 


„  ,  2  2,1/2 

4ir(a  +-  ) 

P 


I(<j>1)  [a2cos  (4>-4>1)  r  p  ]  f ^ (R)  d^1 


.  -co 


i  (41) 


a2sin  (4>-(j)1)  +  P  [4>-4)1-i'4]‘ 


f2(R)  d$J 


>(2.30) 


f  A 


p  = 


2  ir 


The  expression  may  be  written  as 
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\ 


R2  =  2a2  [1  -  cos  (4>  -  4*1)  ]  +  P  (<j>  -  (j)1  +  C)2 


-  lkR  r~ 

ikR  r 

r-  L1  - 


(ik  +  5) 


f2(R)  = 


3  ,  l\ 

A  (lk  +  ^ 


oo 


En<4>,  = 


-irt^rvn  J  1(4,1)<{^l2  +  (^X)4P  J 

-00  ' 


.  ,  if.  ,  C,  —  V 

4lTk  (a  +p  ) 


.  2  I  -lkR 

l'-^e  d(j)l 


R 


(2 .31) 


and,  integrating  by  parts  twice,  the  differentiation  may  be  thrown  over 
to  the  current  as 


En(*,  6)  = 


-it# 

,,2,2  -z.l/2 
4irk  (a  +p  ) 


-ikR 


.U  - 


+  k  I(<j>  )  [a  cost^^  )+p  ] 


d<f> 


(2.32) 


dl 

where  I  and  —  at  +  oo  are  zero. 
d<T 

Equation(2 .31)  is  written  symbolically  as 


11 


oo 

(<t>,  =  J  It*1)  Z«|>-4>1,C)dcj>1 

-00 


(2.33) 


i 

n 
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The  Fourier  Transform  may  easily  be  taken  since  the  integral  is  of  the 
convolution  type.  The  transform,  denoted  by  with  respect  to 


e<B,£)  =  l(B)  Z(z,C) 


(2.34) 


B  is  the  transform  variable  and  b  is  to  be  treated  as  a  parameter., 
-itjff 


If  the  factor 


.  2  2,  1/2 

4 tr(a  +R  ) 


is  ignored  and  z  =  p  tj>  in  the  non-periodic 


part,  Z(4>;£)  may  be  written  as 


„  2  -lkR 

Z(<(>,z,Z)  =  [a  cos  <f>  *•  P  ]  —  - — 


,2  -lkR 
d  e 

~2  R 
dz 


(2  35) 


Note  that  when  a  derivative  with  respect  to  z  is  taken,  4“  is  a  function 
of  z  Now  Z(4>,z,£>)  is  periodic  in  4*  with  period  2 IT.  If  Z(4\z,?=)  is 
expanded  in  a  Fourier  series  with  z  and  as  parameters  the  result  is 


z(.4>,z,&) 


m<j> 

e 


(2.36) 


where 


21T 

Zn(z,^)  =  “  J  Z(<j>,  z,4)  e"ln4>d<j>  (2.37) 

0 


If  z  is  replaced  by  ptj3  in  Equation  (2.36)  the  result  is 


ZC^'pPS)  =  E  Zn(p<J>,t)  ein<^ 
n--co 
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If  the  Fourier  Transform  of  Z  (p4>,£)  with  respect  to  is  Z  (B,4)  the 

n  ’  n 

transform  of  Z(4>,  p4*^)  is 


oo 

Z  (S.C)  =  L  Z  (B+n) 
y  n 

n--oo 


(2.38) 


with 


Z  (8,4) 

n 


(2.39) 


If  Z  (p<}>4)  is  replaced  by  using  Equation  (2.37)  one  obtains,  after  an  inter- 
n 

change  in  the  order  of  integration, 


27 J  oo 

Zn(B,S)  =  (~)  Jf*  e"ln<t>  (ic^1  J  Z(4>1,i4>^)elStf>d4> 


(2.40) 


The  integration  with  respect  to  $  may  now  be  effected  by 


oo 


..  ,  2  2.1/2 

-ik(g  +o  ) 


0/  2  2.1/2 

2(g  +o  ) 


e2B'dg  =  Ko[a(E2-k2)1/2J 


(2.41) 
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using  Campbell  and  Foster  pair  917. 

(, 

Making  the  change  of  variable  a  =  4*  +  — 

P 

1  2  1  -2 
Z ( 4*  ,  a)  =  [a  cos  4  +p  ] 


in  Equation  (2.35)  resul  i.s 


-ikR 

e 

R 


in 
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If 

w=</2  z  (1  -  cos  <f>)1'^(B2-p  k2)1,/2 


then  z  is 


a  ,„2  -  ,2  1/2 

z  =  -  (B  -  p  k  ) 


If  Equation  (2 . 44)  is  substituted  i  nto  integral  (2 . 4(5)  and  integrated,  the 

result  for  Z  (B,C)  is 
n  ’ 

VB’l>  .  *2.  P  §  V'n^lVl1 

P 


(2.45) 


where 


and  all  arguments  are 


1  -2  S2 

B  =  \  (P  -  V 
a  k 


a  ,q2  -  2.  1/2 

—  (.6  -  p  k  ) 


Therefore  from  Equations  (2 . 38)  and  (2 . 45)  Z  (6,^)  is 

i 


Z(B,£)  =  —  E  e 
p  n=-oo 


2  x 


-i (B+n) 


[I  K  (T  )  -f  2BI  K  (T  )  +  I  K  <T  )]  (2.46) 

n-fl  n-*-l  n  nnn  nn-ln 


where 


a  r  to  \2  -  ,  2  •.  1/2 

=  -  [ (S+n)  -  p  k  1 


27 


28 


derived  by  applying  boundary  conditions.  Patton-*-^  obtained  the 
determinantal  equation  for  the  round  wire  version  of  the  helix  by  using 
the  transform  techniques.  The  approximations  used  for  the  round  wire 
version  were  that  the  current  is  concentrated  in  a  filament  at  the  center 
of  the  wire  and  that  the  electric  field  is  zero  on  a  line  which  is  either 
the  smallest  or  largest  distance  from  the  z  axis.  The  transform 
technique  had  not  previously  been  applied  to  the  tape  helix.  The 
approximations  used  were  that  the  current  is  concentrated  in  a  filament 
at  the  center  of  the  tape  and  that  the  average  of  the  electric  field 
over  the  tape  is  zero. 

2.4  Interpretation  of  the  Determinantal  Equation  as  a  Complex-Valued  Equation 

In  addition  to  the  real  valued  solutions  for  the  determinantal 

7 

equation  found  by  Sensiper  it  will  be  shown  in  Chapter  4  that  there 
exist  solutions  where  the  phase  constant,  B,  is  complex-valued  if  the 
determinantal  equation  is  interpreted  as  a  complex-valued  equation.  The 
interpretation  of  the  determinantal  equation  as  a  complex-valued 
equation  is  now  considered. 

If  there  exists  a  solution  for  B  which  is  complex-valued,  then  the 

arguments.  T  of  the  I  K  products  are  complex-valued  since 
m  n  n 

_  r/n  N2  -2, 2,1/2  a 

T  =  [ (B+m)  -  p  k  ]  — 

m  '  - 


The  question  of  which  branch  of  the  square  root  to  choose  for  T  will  be 

m 

discussed  shortly.  Again  consider  the  deteim  inantal  equation.  If  the 
arguments  of  the  IK  products  are  complex-valued,  then  the  IK  products  are 


V 
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complex- valued.  It  is  seen  that  if  8  is  complex-valued  then  ‘I  is  complex- 

valued  and  then,  in  turn,  the  IK  products  are  complex-valued.  Each  side 

of  the  determinantal  equation,  Equation  (2,49)  is  complex- valued,.  It 

remains,  of  course,  to  show  that  there  exist  solutions  for  8  which  are 

complex-valued.  This  will  be  done. 

Which  branch  of  the  square  root  for  ^  should  be  chosen?  If  T 

m  n 

written  as 


T 

n 


>'„2>1/2 


-  T 

nr 


1  T 

ni 


where 


0  =  tan  1 

n 


T 

ni 

T 

nr 


the  question  is  reduced  to  choosing  the  appropriate  the  branch  of  the 
arc  tangent  for  9^„  One  set  of  solutions  of  the  determinantal  equation 
was  found  choosing  the  principle  value  of  the  arc  tangent,  for  each  9  . 

A  second  set  of  solutions  was  found  by  adding  -•  tr  to  the  principle  value 
of  the  arc  tangent  for  0  and  choosing  the  principle  value  of  the  arc 
tangent  for  9^,  n  £  -1.  The  or  -  sign  is  chosen  to  make  larg  T  <  ir. 

A  third  set  of  solutions  was  found  by  adding  r  f  to  the  principle  value 

of  the  arc  tangent  of  both  9  ^  and  9  and  choosing  the  principle  value 

of  the  arc  tangent  for  9  n  ;£  -1,  ~2„  The  f  or  -  sign  is  chosen  to  make 

larg  T  ^ I  <  it,  n  =  -1,  ~2.  For  a  set  of  solutions  each  argument  is  treated 
as  a  continuous  function  of  k.  To  understand  this  continuous  treatment 


t  =  [CB  -l)2  -  (k2  +  8  2)  4-  2i.  8  (8  -1)11/2  - 

-I  r  t  t  v 


consider 


with 


2  B. (B  -1) 

i  r _ 

o  T2  o 
(6  -1)  -(k  +3  ) 
r  x 


where 

i®  , 

T_1  =  I T_1 I  e  ,  k  =  k  p 


and 


B  =  B  +  l  B 
r  i 

Suppose  the  k-B  diagram  has  a  solution  k  =  B^  sin  ^  with  B^  <  0  and  small. 

Figure  2.2a  and  2b  shows  the  k-8  diagram.  The  salient  features  of  the 
2  - 

locus  of  t  as  a  function  of  k  are  shown  in  Figure  2.3.  For  k  >  0 

2 

2  2—2 
and  small  the  real  part  of  T  1  is  positive  since  (B^-l).  >  k  +  B^ 

_  2 

The  imaginary  part  of  T  is  also  positive  since  B^CB^-l)  is  positive. 

2  2 

With  both  the  real  and  imaginary  parts  of  T  1  positive  the  T  ^  as  a 
function  of  k  lies  in  the  first  quadrant.  The  locus  remains  in  the  first 
quadrant  until  a  value  of  k  is  reached  such  that 

2  -2  2 

(B  -1)  =  k  +  B. 

r  l 

This  value  of  k  is  approximately  k  =  1-B^  since  B^  is  small.  As  k  is 
increased  the  locus  enters  the  second  quadrant  to  remain  until 

k  =  sin  41,  i.e.  B  = 

’  r 


1 
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Now  as  k  is  increased  the  locus  of  T  '  enters  the  third  quadrant 

2 

since  8.  <  0  and  6  -1  >  0.  As  it  is  further  increased  <8  -l)2~k  rB 

i  r  r  i 

and  the  locus  leaves  the  third  quadrant. 

Now  for  k  increasing  the  locus  enters  the  fourth  quadrant  since 
2  -2  2 

(8^-1)  >  k  +  8^  and  8X  0.  The  other  arguments  are  treated 

similarly  as  continuous  functions. 

2,5  Summary 

The  same  determinantal  equation  for  the  tape  helix  has  been  derived 
by  two  methods  making  slightly  different  assumptions.  In  the  derivation  by 
matching  boundary  conditions  a  uniform  current  over  the  entire  width  of 
the  tape  was  assumed.  For  the  derivation  by  Fourier  Transforms  a  current, 
was  concentrated  in  a  filament  in  the  center  line  of  the  helix,  and  E 

11 

on  the  tape  was  averaged  to  zero. 

If  the  phase  constant  8  is  complex-valued  then  the  determinantal 
equation  is  complex-valued.  Further,  if  the  phase  constant  is  complex^ 
valued  one  knows  how  the  arguments  of  the  IK  products  can  be  treated  as 
continuous  functions.  The  determinantal  equation  is  now  ready  for  study. 
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3  o  STUDY  OF  TWO  SIMPLIFIED  EQUATIONS 


3.1  Introduction 

The  determinantal  equation,  Equation  (2.26),  is  somewhat  complex  as 
well  as  being  complex  valued  if  the  phase  constant,  S,  is  complex  valued. 
It  will  be  instructive  to  study  two  simplified  equations  so  that  one  will 
gain  an  insight  into  the  behavior  of  the  determinantal  equation.  This 
study  will  be  done  in  the  next  two  sections. 

3.2  The  First  Simplified  Equation 

The  determinantal  equation  with  k  =  k  p  may  be  written  as 


I32  - 


sin2'^ 


00 

E  [I 

n=-oo 


K  .  (T  )  4.  I  K  (T  )  -  2  I  K  (T  )]  D 
nil  n+l  n  n-1  n-1  n  n  n  n  n 


k  ctn 


oo 

2  E 


(3.1) 


n=-oo 


I  K  (T  )  D 
n  n  n  n 


where 


T  =  -  [(B+n)2 
n  - 
P 


k 


2 

] 


1/2 


and 

6 

sin  (B+n)  — - 
D  =  _ 2p 


6 

(Bin)  — 

2p 


The  products  all  have  branch  points  at  the  origin.  To  see  this  consider 
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a  small  argument  approximation  for  the  hyperbolic  Bessel  functions. 
For  | z |  «  1 


I  (z) 
n 


n  ,i  0 


I  (z)  ~  1  K  f'z)  -  In  — 
o  —  o  rj  2 


K  (z) 


n-2 


In 


1  (-)  ( n-1 ) !  of  (1)  ]  nr10 

2  Z  Z  ' 


if 


n  ;£  0 

lim  I  (T  )K  (T  )  -  — —  a  constant 

n  m  n  m  2  n  ' 

T  — »0 
m 

However,  the  first  term  in  the  small  argument  approximation  for  K  (z) 

n 

is  logarithmic  and  has  branch  point  at  the  origin.  This  branch  point 

remains,  although  in  the  limit  as  z— >0  the  term  containing  the 

logarithmic  factor  tends  to  zero.  For  n  =  0  the  lim  I  (T  )K  (t  )  has 

T _ mom 

logarithmic  character.  Consequently,  if  the  argument  of  is  small, 

that  term  will  make  a  substantial  contribution  to  the  determ.inantal 
equation. 

If  one  inspects  T 

2  1 /2 

t  =  —  [B  -l)2-(B.2+k  ) *-1 2B . (B  -1)]  ' 

— 1—  r  i  .ir 

P 


where  8-2  -  3.8  it  is  seen  that  in  the  region  where  B  is  near  unitv 

r  1  r  ’ 

T  will  change  rapidly,,  These  facts,  that  T  changes  rapidly  when  B 
- 1  - 1  r 

is  near  unity  and  that  if  the  argument  of  T  K  will  make  a  substantial 

o  o 

contribution  to  the  determmantal  equation,  suggesi  that  a  simplified 
equation  that  would  behave  like  the  determi nant a  1  equation  near  B  -  1 
would  be  obtained  by  replacing  the  sum  in  the  numerator  of  Equation  (3.1.) 
by 


1  K  (T 
o  o 


r  k  i  f 
2  2 


Inspection  of  the  denominator  of  Equal  ion <3. 1 )  shows  That  the  denominator 
is  slowly  varying  with  6.  This  additional  result  leads  one  to  the 
following  simplified  equation 


I  K  (  T  ,)  ~  T  K  t'T  ) 
go  -1  2  2  -1 


-  a  i 

k  ctn  t 


2  .1  K  <T  ) 
1  1  -1. 


(3, 


where  A  is  a  parameter. 

The  parameter  A  is  picked  to  be  about  equal  to  2A  for  the  values  of 

o 

the  helix  parameter  studied. 

As  the  tape  is  made  more  narrow  the  sum  A  increases. 

o 

The  solutions  for  B  as  a  function  of  k  are  shown  in  Figures  3.1a 
through  3. 10b  for  various  values  of  A.  The  rea.1  and  imaginary  parts  of 
the  phase  constant  are  given  on  individual  figures  with  the  real  part 
as  a  function  of  k  on.  the  figures  wi ' h  suffix  "a"  and  the  imaginary  part 
as  a  function  of  k  on  the  figures  with  suffix  "b”. 
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The  line  k  =  6  sm  is  shown  on  all  the  k  -  B  planes.  The  line 
r  r 

separates  the  k  -  B^_  plane  into  two  regions,  one  the  region  above  the 

line  k  =  sin  ^  where  the  solutions  correspond  to  fast  waves,  i.e. 

waves  which  have  a  phase  constant  corresponding  to  a  phase  velocity- 

greater  than  velocity  of  light,  and  two  the  region  below  the  line  where 

the  solutions  correspond  to  slow  waves,  i.e.,  waves  which,  have  a  phase 

constant  corresponding  to  a  phase  velocity  less  than  the  velocity  of 

light.  The  lines  B  -  1  +  k  are  also  shown  on  the  k  -  B  planes. 

r  r 

Inside  the  "v"  shaped  region  formed  by  the  lines  B^  =  1  *•  k  the  argument 
t  has  non-zero  imaginary  part.  Outside  the  "v"  shaped  region  formed 
by  the  lines  B^  =  1  +  k  the  argument  T  ^  may  have  zero  imaginary  part. 

The  immediate  observation  upon  examining  Figures  3.7a  and  3.7b  is  that 
for  fixed  k  there  are  several  values  of  B.  Some  of  the  solutions  correspond 
to  choosing  7  ^  as  the  argument,  of  rhe  IK  products  and  others  to  choosing 
-  T  The  solutions  are  named  mode  1,  mode  2  and  mode  3  with  the 
determination  of  the  argument  for  mode  1  as  +  T  and  the  determination 
of  the  argument  for  mode  2  and  mode  3  as  -  T  Let.  k  be  the  smallest 
value  of  k  such  that  mode  1  has  a  solution  for  B  with  a  non-zero  imaginary 
part . 

The  mode  1  solution  for  k  <  k  agrees  with  the  solution  shown  by 
7 

Sensiper  for  the  real  valued  determinants!  equation.  The  mode  2 

solution  also  agrees  with  the  solution  found  by  Sensiper. 

The  character  of  the  solution  is  different  depending  on  whether  the 

parameter  A  is  greater  than  A  or  less  than.  A  .  The  mode  1  solution  for 

cl  cl 

k  >  k  is  typical  of  the  solutions  when  A  >  A  where  A  will  be  defined 
cl  c  1  c  1 

later.  The  only  necessary  information  now  is  that  A  used  in  Figures  3.7a 
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and  3.7b  satisfies  A  >  A  ,,  Note  That  real  part  of  the  mode  1  solution 

cl 

is  not  on  the  k  —  8^  sin  ^  line  which  represents  a  solution  corresponding 

to  the  velocity  of  light.  As  k  is  increased  from  k^  it  is  seen  that  8^ 
increases  very  fast  as  a  function  of  k.  The  imaginary  part  of  8  for  mode 
3  is  larger  than  the  corresponding  solution  for  mode  1.  A  wave  with 
this  attenuation  would  decay  rapidly.  Consider  the  results  for  A  =.  1.0 
shown  in  Figures  3.2a  and  3.2b. 

The  solutions  with  zero  imaginary  parts,  mode  2  and  mode  1  with 

k  <  k  ,  are  similar  to  the  corresponding  solutions  where  A  •=  3.0,  This 
cl 

similarity,  however,  is  not  found  for  » he  solutions  with  non-zero  imaginary 

parts.  The  behavior  is  considerably  different  than  when  A  >  A  .  The 

fact  that  there  is  a  considerable  difference  m  behavior  of  the  solutions 

for  A  <  A  ,  and  A  >  A  ,  leads  to  the  obvious  definition  for  A  „ .  The 
cl  cl  cl 

value  A  is  that  value  of  the  parameter  A  in  Equal i on (3 . 2)  that  separates 

the  values  of  A  into  two  sets,  one  set,  A  >  A  ,  such  that  the  solutions 

have  the  character  similar  to  that  shown  in  Figures  3.7a  and  3.7b  and  the 

other  set,  A  <  A  ,  such  that  the  solutions  have  the  character  similar  to 
cl 

that  shown  in  Figures  3.2a  and  3,2b.  The  outstanding  difference  m  the 

behavior  for  A  <  A  and  A  >  A^  may  be  observed  by  inspecting  the  imaginary 

part  of  the  phase  constant  corresponding  to  mode  1.  If  A  >  A  then 

cl 

8.  >0  for  all  k  >  k  .  If  A  <  A  ,  then  there  exist  regions  where  there 
l  cl  cl 

exists  no  mode  1  solution  for  k  >  k  ,  .  The  value  A  =  A  „  corresponds  to 

cl  cl 

a  value  2A^  which,  in  turn,  corresponds  to  a  tape  width  5  that  is  larger 
than  tape  widths  consistent  with  the  narrow  tape  approximation.  Therefore, 
no  solutions  of  the  determinant al  equation  in  the  region  8^  near  unity  will 
correspond  to  the  solutions  given  in  Figures  3.2a  and  3.2b.  They  are  given 


1 


I 

to  assist  in  the  understanding  of  this  kind  of  equation,,  The  strangest  * 

property  of  the  solutions  in  Figures  3.2a  and  3.2b  is  that  in  order  to 

have  a  continuous  solution,  B  as  a  function  of  k,  it  is  necessary  to  [ 

change  branches  of  the  argument  at  k  =  .273  and  k  =  .358.  These 

values  of  k  correspond  to  values  of  8  =  0.  | 

The  argument  for  the  mode  1  solution  corresponding  to  k  =  .372  is  j 

pe ^  2  where  p  is  real.  The  argument  for  the  mode  2  solution  corresponding 
to  it  =  .273+  is  pe^^^2  ^  .  The  arguments  are  discontinuous  by  the  | 

amount  ir  but  the  value  of  the  right  hand  side  of  Equation(3.2}is  the  same 
for  both  arguments.  The  reason  for  this  equality  is  not  at  first  obvious.  | 

Now 


lo(ip)Ko(ip)+I2(ip)K2(ip)-2I1(ip)K1(ip) 


is  the  complex  conjugate  of 

Io(-ip)Ko(-ip)+I2  (-ip)K2(-ip)  -2I1(-ip) K^-ip) 


for  p  real.  In  addition,  the  sum  of  the  imaginary  parts  is  zero  for 

the  value  k  =  .273.  This  may  be  seen  easily  from  the  fact  that  for 

k  =  .273,  B^  =  0,  and  therefore,  the  left  hand  side  of  Equation  (3. 2)has 

zero  imaginary  part.  The  right  hand  side  also  has  zero  imaginary  part. 

The  right  hand  sides  for  k  =  .273  and  .273'  are  complex  conjugates  and  have 

zero  imaginary  parts,  and  therefore  the  right  hand  sides  are  equal.  A 

similar  situation  occurs  for  k  =  .358.  Figures  3.6a  and  3.6b  show  the 

results  for  a  value  of  A  near  A  ,,  A  =  1.5. 

cl- 
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Note  that  the  difference  in  behavior  of  the  solution  for  A  =  1.0  and 

A  =  1.5.  Both  A  =  1.0  and  A  =  1.5  are  less  than  A  ,  since  there  exists 

cl 

a  region  of  k  >  where  no  solui ion  tor  mode  1  exists.  The  real  parts 

of  the  mode  1  solution  are  connected  by  a  continuous  solution  of  the 

other  determination  of  the  square  root  compared  to  the  mode  1  solution 

for  A  =  1.5.  However,  for  A  =  1.0  the  real  parrs  of  rhe  mode  1  solution 

are  not  connected  by  a  continuous  solution.  If  Figures  3„2a  through  3.6b 

are  studied,  it  is  seen  that  the  change  from  connecting  the  real  parts  of 

the  mode  1  solution  by  a  continuous  solution  ’0  no’,  connecting  the  real 

parts  takes  place  in  a  continuous  manner.  Thtre  is  no  quesiion  of 

discontinuity  of  the  arguments  since  all  the  solutions  hare  the  same 

determination  for  the  square  roots.  It  is  reitera'ed  'har  this 

particular  behavior  will  not  be  found  in  the  determinant al  equation 

since  A  <  A  ,  corresponds  to  tapes  which  a.re  too  wide  for  the  narrow 
cl 

tape  approximation. 

Figure  3.1a  and  3.1b  show  the  results  for  A  =  .9  and  it  is  observed 
that  the  behavior  of  the  solutions  as  a  function  of  A  is  in  a  continuous 
manner . 

Consider  some  of  the  general  properties  of  the  solutions  for  A  >  A  .. 

cl 

(Inspect  Figures  3.7a  through  3.10b)  Note  that  k  .  the  lowest  value  of  k 
for  which  mode  1  has  a  non-zero  imaginary  part,  becomes  greater  as  A  is 
increased.  Also  for  k  >  k  the  real  part  of  S  follows  more  closely  to  the 
k  =  sin  ']J  line  and  follows  for  a  greater  range  of  k  values.  For  the 
imaginary  part  the  local  maximum  near  k  =  .2  and  the  local  minimum  near 
k  =  .25  decrease  as  A  increases.  The  values  of  for  ,4  <  k  <  .5 
increase  as  A  increases.  For  mode  2,  the  imaginary  part  increases  as  A 
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Figure  3.2a, 


The  k  -  (3  diagram  for  simplified  equation  one; 
A  =  1.0,  r4  ~  12.6°. 


Figure  3.2b„  The  k  -  B  diagram  for  simplified  equation  one; 
A  =  1,0,  ^  12.6°. 
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Figure  3.3a.  The  k  -  B  diagram  for  simplified  equation  one 
expanded  near  k  =  .25;  A  =  1.0,  =  12.6°. 


0  .02  .04  ,06  .08  .10  .12  .14 

/S'- 


Figure  3.3b.  The  k  -  B^  diagram  for  simplified  equation  one 
expanded  near  k  =  .25;  A  =  1.0,  ^  =  12.6  . 


MODE  3. 


-  T 


Figure  3.4b*  The  k  -  B  diagram  for  simplified  equation  one 
expanded  near  k  -  .25;  A  =?  1,05,  ^  =  12.6°, 
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Figure  3.6a.  The  k  -  B  diagram  for  simplified  equation  one, 
A  -  1.5^  r'F  -  12 o 6 
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Figure  3.6b.  The  k  -  B.  diagram  for  simplified  equation  one., 
A  =  1.5,  =  12.6°. 
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Figure  3.9b.  The  k  -  B  diagram  for  simplified  equation  one 
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increases . 

One  could  study  this  simplified  equation  in  more  detail,  e.g.,  as  a 
function  of  pitch  angle  41.  This  kind  of  study  will  be  delayed  until  the 
complete  determinantal  equation  is  studied.  At  that  time  the  physical 
significance  of  the  increasing  or  decreasing  of  solutions  will  be 
discussed. 

The  reason  for  studying  the  simplified  equation  was  to  show  the  exist¬ 
ence  of  the  various  modal  solutions  of  the  determinantal  equation  and 
show  how  the  arguments  of  the  IK  products  could  be  treated  so  that  the 
solutions  would  be  continuous  functions  of  k. 

3.3  The  Second  Simplified  Equation 

It  will  be  instructive  to  study  another  simplified  equation.  This 

equation  will  approximate  the  behavior  of  the  determinantal  equation  in 

the  region  8  near  2.  Since  for  8  ~  2, 

r  r  ~  ’ 

2 

T  0  -  f (B  -2) 2  -  (k  +B.2)+2iB  (8  -2)]1/2ctn  t 

_ ■  2  r  i  x  r 


changes  rapidly,  only  those  terms  in  the  numerator  of  the  determinantal 
equation.  Equation  0. 1) ,  will  be  used  for  the  simplified  equation.  Also, 
as  stated  in  the  last  section,  the  denominator,  2Aq,  is  a  slowly  varying 
function  of  k  and  8.  Again  2A^  is  replaced  by  a  constant,  A.  Simplified 
equation  two  is 


B2 - — 


sin 


WV  +  WV  -2  WV 


k  ctn  ^ 


(3.3) 
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The  solutions  are  shown  in  the  following  figures. 

The  line  k  =  B  sin  4*  is  shown  on  the  k  -  B  planes.  The  lines 
r  r 

B  =  2  v  k  are  shown  on  the  k  -  B  planes.  Inside  the  *vM  shaped 
r  -  r 

region  formed  by  the  lines  B^  =  2  -i-  k  the  argument  "T  has  non-zero 
imaginary  part.  Outside  the  "v"  shaped  region  formed  by  the  lines 
B  =  2  4-  k  the  argument  t  ^  may  have  zero  imaginary  part. 

If  one  examines  Figures  3.11  and  3.12  the  immediate  observation  is 
that  for  a  given  k  there  exists  more  than  one  solution  for  B.  Tne  different 
modes  do  not  ail  have  the  same  determination  of  the  square  root  for 
agrument  of  the  IK  products.  Modes  2  and  5  have  the  positive  deter¬ 
mination  for  the  argument  of  the  TK  products,  whereas  modes  4  and  6  have 
the  negative  determination. 

Let  k  „  be  the  lowest  value  of  k  such  that  the  solution  for  mode 
c2 

2  will  have  a  non-zero  imaginary  part  .  The  mode  2  solution  for 

B  <  2-k  and  k  <  k  is  similar  to  the  solution  found  by  Sensiper7 
r  c2 

The  mode  4  solution  also  corresponds  to  Sensiper’s  results.  The  mode  2 

solution  for  k  >  k  has  an  imaginary  part  which  is  similar  to  the  mode 

one  solution  to  the  simplified  equation  of  the  previous  section. 

Consider  the  results  shown  in  Figure  3.16  which  is  an  expanded  view 

of  a  section  of  Figure  3.11a„  Note  that  to  have  the  real  part  of  the 

phase  constant  continuous  across  the  line  B  =  2-k.  for  k  <  k  ,  it  is 

r  c2 

necessary  to  change  The  determination  of  the  square  root.  This  difference 
in  the  determination  of  the  arguments  causes  a  discontinuity  in  the 
arguments  by  the  amount  IT.  However  the  value  of.  the  right  hand  side 
of  Equation  0*3)  is  the  same  for  either  determination  of  the  T 
argument  at  B^  -  2-k.  This  equality  is  true  because  when  B^  =■  2-k 
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and  B  =0,  T  „  =  0  and  further 

l  ’  -2 


lim  In(T-2)  W  =  aS  ’  "  *  ° 

T-2  *  ° 


which  is  independent  of  which  branch  was  chosen. 

The  results  shown  in  Figure  3.18  which  is  an  expanded  view  of  a 

section  of  Figure  3.11a  are  to  be  studied.  Again  the  same  type  of 

discontinuity  occurs.  The  line  Bf  =  2  +  k  separates  inodes  4  and  5  which 

have  a  continuous  real  part.  The  argument  has  a  discontinuity  equal  to  v 

for  modes  4  and  5  at  B  =  2  +  k.  The  right  hand  side  of  Equation  3.3 

r 

is  the  same  for  either  determination  of  the  argument.  The  reason  for  this 

equality  is  similar  to  that  given  above. 

In  Figure  3.21a  the  behavior  of  the  solutions  for  modes  2  and  5  should 

be  observed.  This  behavior  is  typical  for  values  of  the  parameter  A  less 

than  A  .  where  A  will  be  defined  later.  In  Figure  3.19  the  behavior 
c2'  c2 

of  the  solutions  for  modes  2  and  5  for  A  =  10. should  be  studied.  Note 

the  behavior  of  the  solutions  is  different  for  A  =  10.0  than  for  A  =  7.0. 

This  behavior  for  A  =  10.0  is  typical  for  values  of  the  parameter  A  >  A  . 

The  number  A  is  now  easily  defined.  The  number  A  _  divides  the  set  A 

c2  c2 

into  two  classes,  one,  A  <  A  ,  where  the  mode  2  solutions  corresponding 
to  these  values  of  A  are  not  continuous  with  the  solutions  for  mode  4 
and  the  other,  A  >  A  ,  where  the  mode  2  solutions  corresponding  to 

C  j£ 

these  values  of  A  are  continuous  with  the  solutions  for  mode  4. 

If  one  examines  Figures  3.22a,  3.22b,  3.23a,  3.23b, .3.24a,  and  3.24b 
it  is  seen  how  these  solutions  take  on  the  different  character  as  the 
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Figure  3.15a.  The  k  -  B  diagram  for  simplified  equation  two; 
A  =  1.0,  r+  —  12.6°. 


Figure  3.15b.  The  k  -  B.  diagram  for  simplified  equation  two; 
A  =  1.0,  14J  =  12.6°. 


Figure  3.16.  The  k  -  B  diagram  for  simplified  equation  two  near  k  ;  A  =  1.0  and  5.0 


Figure  3.17.  The  k  -  8r  diagram  for  simplified  equatio 
two  near  k  A  =  10.0  and  25.0,  41  -  10° 


Figure  3.19.  The  k  -  Si  diagram  for  simplified  equation  two  near  k  •  A  =  10.0  and  25 


Figure  3.20a.  The  k  -  B  diagram  for  simplified  equation  two  near  k 
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Figure  3.24a.  The  k  -  8  diagram  for  simplified  equation  two  near  k 
A  =  8.0,  ^  =  10°. 


Figure  3.24b. 


0i 

The  k  -  B.  diagram  for  simplified  equation  two  near 
A  =  8.0,  ^  =  10°. 
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value  of  A  is  increased  through  the  value  A  .  There  is  no  question  of 

c2 

discontinuity  in  the  argument  here  since  all  arguments  are  of  the  same 
determination.  Note  that  the  character  of  the  solutions  changes  in  a 
continuous  manner  as  a  function  of  the  parameter  A.  This  type  of 
behavior,  as  well  as  the  behavior  exhibited  at  the  B  =  2  +  k  lines, 
will  be  seen  in  the  determinantal  equation. 

3.4  Summary 

Two  equations  which  are  relatively  simple  compared  to  the  determinantal 
equation  and  which,  furthermore,  exhibited  some  of  the  character  of  the 
determinantal  equation  were  studied  in  this  chapter.  It  was  shown  how 
discontinuities  in  the  argument  of  IK  product  could  result  in  solutions 
which  would  be  continuous.  If  the  discontinuities  of  the  arguments  are 
not  introduced  there  is  no  continuation  of  the  solutions.  Further  it 
was  seen  that  the  solutions  for  modes  2  and  4  changed  character  as  a 
function  of  the  parameter  A. 

The  parameter  A  was  chosen  to  be  approximately  equal  to  2Ao*  As  the 
tape  width  5  is  made  more  narrow  A^  increases.  Consequently,  the  larger 
values  of  the  parameter  A  correspond  to  The  narrower  tape  widths. 

Solutions  for  S  which  are  complex  valued  have  been  found  for  both 
simplified  equations.  The  use  of  complex  roots  will  be  discussed  later. 

I 

The  behavior  of  the  solution  of  simplified  equation  one,  Equation  (3.2), 

for  A  =  1.0  is  reminiscent  of  Pierce's  coupled  mode  theory  for  the  helix 

25 

as  discussed  by  Watkins  .  However,  as  was  pointed  out  in  the  discussion 

of  equation  one,  the  values  of  the  parameter  A  <  A^  correspond  to  tapes 

which  are  not  narrow,  i.e.  too  wide  for  the  determinantal  equation  to  be 

valid.  The  fact  that  the  solutions  for  A  >  A  , ,  which  correspond  to 

cl 
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narrow  tapes,  do  not  exhibit  the  character  needed  for  Pierce's  coupled 

mode  theory  inhibits  the  use  of  Pierce's  simple  theory. 

Simplified  equations  one  and  two,  Equation  (3.2.)  and  (3.3),  are 

related  to  the  determinantal  equation  for  the  sheath  helix  for  sheath 

helix  modes  -1  and  -2  respectively.  The  complex-valued  solutions  for 

the  sheath  helix  and  a  discussion  of  their  relation  to  the  solutions 

29 

for  the  tape  helix  are  given  by  Klock  and  Mittra 

With  the  background  obtained  from  studying  the  solutions  of  the 
two  simplified  equations  one  is  now  ready  for  the  determinantal  equation. 
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4.  SOLUTION  OF  THE  DETERMINANTAL  EQUATION 


4.1  Introduction 


The  determinantal  equation 


sin2^ 


k  ctn  4* 


A  +A  -2A 
+  1  -1  o 

2A 


(4.1) 


oo 

E 

n=-oo 


I  .(T  ')  K  (T  )  D 
n-»-j  n  n+j  n  n 


2  -2  1/2  a 

Tn  =  [(6  +  n)2  -  k  ]  § 

P 


sin  (Bt-n)  - 


D  = 
n 


2p 


,  k  =  k  p 


(34-n)  ^ 


derived  in  Chapter  2  which  is  a  complex-valued  equation  if  B  is  complex- 
valued.,  was  solved  with  the  aid  of  a  digital  computer. 

The  determinantal  equation  can  be  written  such  that  the  complex-valued 
phase  constant;  B;  is  a  function  only  of  the  parameters,  normalized  tape 
width;  pitch  angle,  and  normalized  frequency.  Symbolically,  this  is 
written  as  B(^;4*;  k)  .  The  tape  width  is  normalized  with  respect  to  the 

r  6  a  i 

radius;  i.e,  0  =  — .  The  ratio  —  is  replaced  by  its  equal,  ctn  V;  so  that 

p 

Equation  (4.1)  reads 


B2 - — 


Cd 

k  ctn  4* 


A  +A  -2A 
4-1  —  1  o 

2A 


(4.2) 
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A 

J 


CO 


I  .  (T  ) 
n+j  n 


K  .  f  T  ) 
n+j  n 


D 

n 


2  1/2 

t  =  [(B+n)  -  k  ]  ctn  41 

n 

sm  [  (B+n)  —  ctn  41] 

D  = - - - 

n  (B+n)  _  ctn  4' 

2 

6  . 

Define  a  -  —  ctn  ^ 

_  ks 

The  normalized  frequency  is  k  =  ------qj  where  k  is  the  free  space  wave 

number.  By  normalizing  with  respect  to  the  radius  one  is  able  to 

reduce  the  study  of  the  determinantal  equation  to  the  solution  for  the 

complex  valued  phase  constant  as  a  function  of  only  normalized  tape  width, 

normalized  frequency  and  pitch  angle. 

All  of  the  coefficients  in  Equation  (4.2)  are  real.  Therefore,  if 

B  =  B  +iB  is  a  solution,  then  B  =  B  ~iS  is  a  solution.  If  the  solution 
r  i  r  l 

for  the  phase  constant  is  B  =  B  -lB  ,  then  the  <j>  dependent  fields  are 

~i6  4>  -B  <j>  r  1 

r  l 

of  the  form  e  e 

The  real  part  of  the  phase  constant,  B^,  is  a  measure  of  the  phase 
delay  as  the  wave  travels  in  the  +<j>  direction  and  has  units  of  radians 
per  radian.  The  imaginary  part  of  the  phase  constant,  B^  is  a  measure 
of  the  attenuation  as  the  wave  travels  in  the  +4>  direction  and  has  units 
of  nepers  per  radian. 

The  line  k  =  8  sin  4*  is  shown  on  all  the  curves  for  8..  This  line 
r  l 

separates  the  k  -  B^  diagram  into  two  regions: 


1)  Points  above  the  line  represent  waves  which  are  traveling  down 
the  tape  with  a  phase  velocity  greater  itian  i.he  velocity  of  .light 
i.e.  fast  waves  and 

2)  Points  below  the  line  represent  waves  which,  are  traveling  down 
the  tape  with  a  phase  velocity  less  than  the  velocity  of  light, 
i.e.  slow  waves. 

The  lines  6  =  1  +  k  and  S  •=•  2  +  k  are  shown.  Between  the  lines 

r  r 

8  =  1  +  k  and  between  the  lines  8  =  2  t  k  the  arguments  T  and  T 

r  -  i  -  -1-2 

2  —  ^  1/2  a 

respectively  have  non-zero  real  parts.  Recall  T  =  [ (B+nl  -k  ]  '  —  . 

P 

4.2  Solutions 

4,2.1  The  k  -  8  diagram-  41  =  1 0°,  &  =  .035 

Figures  4.1a,  b,  c,  d  and  e  is  the  k  -  8  diagram  for  pitch  angle  10 
and  5  =  .035.  The  immediate  observation  is  that  for  a  value  of  k  there 
exist  more  than  one  solution  for  8,  The  various  solutions  or  modes  each 
correspond  to  a  different  determination  of  a  square  root  argument  for 
the  IK  products. 

The  determination  of  the  square  roots  was  discussed  in  Section  2.4 
Note  that  mode  1  has  positive  determination  for  both  T  ^  and  T  All 
the  other  modes  have  the  negative  determination  for  T  .  Modes  2,3, 
and  5  have  the  positive  determination  and  modes  4  and  6  have  the 
negative  determination  for  T 

In  the  regions  8^  ~  1  and  8^  ~  2  the  solution  behaves  similarly  to 
the  solution  of  the  simplified  equations  studied  in  Sections  4.2  and  4.3 
respectively. 

The  constants  k  ,  and  k  are  defined  as  in  Section  4.2  and  4.3 
cl  c2 


respectively. 


Figure  4!>]a<>  The  complete  k  “  B  diagram  for  the  tape  helix: 


74 


I 
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Figure  4.1b.  The  k  -  B.  diagram  for  the  tape 
helix:  Modes  2  and  3;  4*  •=  10,  < 


Figure  4„ic„  The  k  —  S  diagram  for  the  tape  helix"  Mode 


Figure  4.1e0  The  k  -  0  diagram  for  the  tape  hell  a  Modes  4  and  5;  r  r  10, 
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The  solutions  for  mode  1,  k  <  k  , .  mode  2.  k  <  k  ;  and  mode  4. 

cl  c2 

k  <  k  •  correspond  and  agree  with  the  results  obtained  by  Sensipe.r. 
c4 

All  other  solutions  have  non-zero  imaginary  parts, 

4.2.2  Mode  1  Solution.  ^  -  12 .6°,  a  =  10~3°8;  .0001,  .1,  . 4 , 

The  solutions  for  mode  1  as  a  function  of  6  are  shown  in  Figures  4„2a 

and  b„  As  the  tape  is  made  more  narrow  the  real  part  of  the  phase  constant 

tends  toward  the  k  =  sin  41  line  as  would  be  predicted  for  an  infinitesimally 

narrow  tape  or  filament.  The  maximum  value  of  the  imaginary  part  of  the 

phase  constant  increases  as  the  tape  is  made  wider.  The  imaginary  part 

of  the  phase  constant  increases  rapidly  as  k  is  increased  from  k  .  The 

cl 

behavior  of  the  solution  changes  quite  fast  for  a  small  change  in  normalized 
frequency. 

Expanded  views  for  a  range  of  k  near  k  are  shown  in  Figures  4,3a 

and  b.  As  the  tape  is  made  more  narrow,  k  ,  increases  with  k  approaching 
_  'cl  c 1  b 

sin  ^  fi  , 

- : — nr  as  0  tends  toward  zero. 

1+sin  V 

4.2.3  Mode  1  Solution,  &  =  ,035,  various  pitch  angles 

The  solutions  for  mode  1  for  various  pitch  angles  for  fixed  tape  width, 

6  =  ,035,  are  shown  in  Figures  4,4a  and  b.  As  the.  pitch  angle  is  increased 

the  S  locus  deviates  less  from  the  k  =  B  sin1!1  line.  This  result  seems 
r  r 

reasonable  physically  since  as  the  pitch  angle  is  increased  a  section  of 
unit  length  of  the  helix  is  becoming  more  nearly  a  straight  wire.  Also 
note  that  as  pitch  angle  is  increased  the  imaginary  part  of  the  phase 
constant  decreases, 

4.2.4  Mode  3  Solution:  ^  =  12,6°,  a  =  10  308 ,  .0001,  ,1,  ,4 
Figures  4.6a  and  b  shows  the  mode  3  solutions  as  a  function  of  tape 

width.  Note  that  the  imaginary  part  of  the  phase  constant  is  now  to  a 


Figure  4.2b.  The  k  -8.  diagram  for  the  tape  helix:  Mode 


Figure  4.3a.  The  k  ~  8  diagram  for  the  tape  helix  near 
Mode  1;  r  *|»  =  12.6°,  a  -  10~308  .0001, 


329 


Figure  4.4b,  The  k  -  0.  diagram  for  the  tape  helix:  Mode 
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Figure  4.5b.  The  k  -  B.  diagram  for  the  tape  helix:  Mode  1;  °  =  .035 


Figure  4,6a0  The  k  6r  diagram  for  the  tape  helix  Mode 
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Figure  4,6b.  The  k  -  B  diagram  for  the  tape  helix.  Mode 
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different  scale  since  G^  is  larger  than  it.  is  for  any  other  mode. 

4.2.5  Modes  2  and  6 .  4*  -  12.6,  a  ^  10  ,  .0001,  .1,  .4 

The  solution  for  modes  2  and  6  are  shown  in  Figures  4.7a  and  b.  Figures 

4.8a  and  b  show  an  expanded  view  for  a  range  of  k  values  near  k  .  The 

behavior  for  mode  2  above  k  is  similar  to  the  behavior  for  mode  1  above 

c2 

k  ..  In  contrast  to  mode  1  behavior,  the  behavior  of  the  solution  with 
cl 

zero  imaginary  part  does  not  approach  the  line  B  =  2  -  k  asymptotically 
but  terminates  at  the  line.  To  find  a  continuous  solution,  the  determination 
of  T  must  be  changed.  This  type  of  behavior  was  exhibited  by  simplified 
equation  two,  Equation  (3.3). 

4.2.6  Modes  4  and  5 :  ^  =  10°,  a  =.  ,10~3°8,  .0001,  .4,  ^  =  12.6°,  a  =,1 

Modes  4  and  5  exhibit,  behavior  as  shown  in  Figures  4„9a  and.  b.  The 

solution  for  mode  4  terminates  on  the  line  B  =  2->-k.  To  have  a  continuous 

r 

solution  the  determination  of  T  2  must  be  changed  as  m  simplified 
equation  two,  Equation  (3.3.1).  Figures  4.10a  and  b  give  the  results  for 
modes  4  and  5  for  different  parameters.  The  solution  exhibits  the  same 
character  as  for  the  previous  parameters. 

4.2.7  Modes  2,4,  and  5r  ^  -  10°,  Br  ~  2.4,  a  =  10~'10, 10~15,  lo"17, 


In  Figure  4.11a  through  Figure  4.15b  are  shown  the  change  of 
character  of  the  solution  that,  modes  2,  4,  and  5  exhibit  as  the  tape  width 
is  changed.  If  the  tape  is  made  sufficiently  narrow,  mode  2  terminates 
on  the  line  B  =  2  t  k  and  connects  in  a  continuous  manner  to  mode  4.  This 
is  exactly  the  same  behavior  as  was  shown  by  simplified  equation  two, 
Equation  (b  >3). 
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Figure  4„7b,  The  k  8,  diagram  for  the  tape  helix  Modes  2  and  6; 

i  o  n  o 


Figure  4.8a.  The  k  -  Br  diagram  for  the  tape  helix  near  k  :  Modes  2  and  6 
^  -  12 , 6°  ,  a  -  10~308  ,0001,  1,  .4, 


Figure  4.9a.  The  k  -  Br  diagram  for  the  tape  helix  Modes  4  and 


Figure  4 ,9b  The  k  -  Bx  diagram  for  the  tape  helix  Modes  4  and 
=  12  6°,  a  =  ,1 


Figure  4.10b,  The  k  -  diagram  for  the  tape  helix  Modes  4  and  5; 
4*  =  10°,  a  =  10-308,  0001,  4 


Modes  2 


Figure  4.13a. 


The 

4, 


k  -  B  diagram  for  the  tape  helix  near  k  • 
and  5;r  +  =  10°,  a  =  lO-”.  °4 


Figure  4.14a.  The  k  -  S  diagram  for  the  tape  helix  near  k  :  Modes  2, 
4,  and  5;^  =  10°,  a  =  KT20.  C4 


Figure  4.14b. 


The  k  -  B.  diagram  for  the  tape  helix  near  k  •  Modes  2. 
4,  and  5;^  =  10°,  a  =  10“20.  c4 
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4.2.8  Solution  when  the  Tape  Width  is  Infinitesimally  Narrow, 
Mittra  ^  showed  that  as  the  tape  width  is  made  infinitesimally 
narrow  the  asymptotic  solution  for  the  determmantal  equation,  Equation 
(2.40)  is  given  by 


3  = 


k 

sin  ^ 


(4.3) 


and 


S  =  1  +  k 


A  graph  of  the  above  solution  is  shown  in  Figure  4.16. 

To  show  that  the  solution  given  in  Equation  (1.3)  is  the  asymptotic 
solution  consider  the  determmantal  equation. 


B 


2 


A  4- A  -2A 
-i-1  -1  o 

2A 

o 


(4.4) 


Recall 


oo 

A.  =  E  I  .K  .(T  )  D 
J  „  „  n+J  n+J  n  n 

n=~oo 


For  an  infinitesimally  narrow  tape  =  1.  Now  for  n  ■£  0 


I  .K  .(T  )  ~ 
n+j  n+j  n  ~ 


JL 

2 


[  (  n+j  ) 


1 

V  2 11/'2 


n 


as  was  shown  by  Sensiper 


7 


If  T  0;  the  numerator  of  the  right  hand 
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side  of  Equation  (4.4)  is  finite.  As  t  ^  approaches  zero  the  numerator  of 

the  right  hand  side  becomes  finite  because  the  I  K  (T  )  term  has  a 

o  o  -1 

logarithmic  singularity  (see  Section  3.2).  The  denominator  of  the  right 
hand  side  of  Equation  (4.4)  for  |n|»l  behaves  like  the  harmonic  series 
and  consequently  the  denominator  of  the  right  hand  side  is  finite.  The 
value  of  the  right  hand  side  of  the  determinantal  equation,  Equation  (4.4), 
is  zero  if  T  ^  0  and  indeterminant  if  T  =  o.  The  determinantal  equation 
then  has  the  asymptotic  solution 


3  = 


k 

sin  41 


and  also  has  the  solution 


T_1  =  0,  i.e.  6  =  1  +  k 

The  solution  obtained  when  6  =  10  a  very  narrow  tape,  agrees  with 

the  asymptotic  solution.  As  the  tape  is  made  wider,  but  still  very  narrow 

k  -30 

e.g.  corresponding  to  0  =  10  ,  the  solution  for  the  phase  constant 

becomes  complex  for  k  >  k  ,  and  the  real  part  of  the  phase  constant 

cl 

deviates  from  the  3  -  — — rr  line  which  is  the  asympte-tic-  solution, 

r  sin  HJ 

For  any  tape  corresponding  to  real  dimensions,  even  very  narrow  tapes,  the 
solution  deviates  from  the  asymptotic  solution. 

4.3  Choosing  the  Predominant  Modes,, 

In  this  section  it  will  be  shown  that  modes  1,  2  and  4  are  permitted 
by  the  interpretation  of  leaky  mode  theory,  whereas  modes  3,  5  and  6  are 
not  physically  admissable. 
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If  only  the  r  and  z  variation  is  considered,  the  fields  outside  'he 
helix  for  r  large  may  be  written  as 


^  (r ,  z) 


■r  * 

oo  n  a 

L  a  e 
n 

n=-oo 


(4.5) 


“■ x 

because  I  K  (x)  behaves  like  e  for  large  x. 
mm 

dependence  is 


Now  if  8  =  8  ~i.B  ,  the  z 
r  v 


which  corresponds  to  a  wave  traveling  in  the  ~z  direction  and  decaying 

in  the  +z  direction.  For  large  r  the  fields  are  given  by  Equation  (4.5). 

The  only  terms  which  could  cause  the  fields  Ti  behave  improperly  are 

the  terms  corresponding  to  n  -  -1,  -2  and  -3.  the  argument  7  7  or 

7  which  could  cause  the  mode  to  behave  improper  lv  is  shown  tor  each 

mode  in  Figure  4.17  along  with  the  r  dependence  for  'he  field.  In  the 

r  dependence  shown  in  Figure  4.17  the  real  members  a  and  b  are  positive. 

ar  >*<  i 

Consider  the  mode  1  solution.  For  8  <1  the  :r  dependence  is  e  e 

r 

which  corresponds  to  an  outwardly  traveling  wave  which  is  attenuated. 

Since  <  1  the  radiation  corresponding  to  the  phase  constant  is  in  the 

backward  direction,  and  decay  in  the  transverse  direction  is  proper  fpr 

18 

the  leaky  wave  as  was  shown  by  Marcuvitz  and  later  discussed  by 

Oliner  ^  .  For  >  1  the  ;r  dependence  is  edTe  1'5r  which  corresponds 

to  an  outwardly  traveling  wave  which  is  growing  in  amplitude.  Since 
8^  >  1  the  radiation  corresponding  to  the  phase  constant  is  in  the 
forward  direction,  and  exponential  growth  in  the  transverse  direction  is 


Figure  4.17a.  Determination  of  T  for  Modes  1  and  2. 


MODE  3 


MODE  4 
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proper  for  the  leaky  wave  as  was  shown  by  Marcuvitz.  To  understand  the 
growth  in  the  transverse  direction  consiaer  a  slotted  waveguide  shown 
schematically  in  Figure  4.18,  The  guide  is  fed  at  the  origin,  and  radiation 
takes  place  in  the  forward  direction.  As  z  is  increased  less  radiation 
occurs  since  the  fields  are  decaying  exponentially.  The  spacing  of  the 
lines  in  Figure  4.18  is  inversely  proportional  to  the  power  density. 

Note  that  at  z  =  z'  the  fields  increase  as  x  is  increased  from  zero  to  a 
value  of  x  related  to  the  direction  of  radiation. 

The  mode  2  solution  has  exactly  t.he  same  behavior  for  T  as  the  mode 
1  solution  has  for  T  ^  and  consequently  mode  2  has  solutions  which 
correspond  to  waves  which  are  leaky  waves.  The  mode  4  solution  corresponds 
to  waves  which  are  leaky  waves  since  the  mode  4  solution  behaves  for  T  ^ 
the  same  as  mode  1  and  mode  2  solutions  behave  for  T  and  T  respectively. 

However,  this  behavior  of  the  fields  is  not  the  case  for  solutions 
corresponding  to  modes  3 ,  5,  and  6.  Each  mode  has  a  transverse  behavior 
which  is  characterized  by  an  inwardly  traveling  wave.  If  there  is  an 
inwardly  traveling  wave  the  field  must  increase  in  the  z  direction.  The 
fields  do  not  increase  since  the  z  dependence  is 

B  B 

r  l 

-l  —  z  ~i  ■-=-  z 

e  5  e  P  B  >  0 

i. 

This  contradiction  between  the  wave  traveling  inwardly  and  the  fields 
decaying  in  the  +z  direction  make  the  waves  corresponding  to  the 
solutions  of  modes  3^  5,  and  6  physically  inadmissable  for  leaky  waves. 

If  B  =  B  +iB.  is  chosen  as  the  solution  to  the  determinantal 
r  i 

equation^  all  of  the  arguments  are  the  conjugates  of  the  arguments  when 
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B  =  3  ~xB  is  chosen  as  the  solution  to  the  deternunant.al  equation, 
r  1 

The  property  that  t  is  its  conjugate  when  B  is  its  conjugate  makes  the 
discussion  of  admissable  solutions  when  B  =  B  --1.B  similar  to  the  above 

t  i. 

discussion  when  B  =  B  -iB  and  with  the  same  results, 

r  l 

1)  Modes  1,  2,  and  4  are  physically  admissable  for  leaky  modes  and 

2)  Modes  2,  5  and  6  are  not  physically  admissable  for  leaky  modes. 

If  -B  is  chosen  as  the  solution  to  the  det erminantal  equation,  then 

attention  must  be  focused  on  the  terms  containing  t  t  ^  and  T  ,  Here, 

as  before,  the  discussion  of  the  admissable  solution  is  similar  to  that 

discussion  above  and  with,  or  course,  the  identical  results.  The  normalized 

frequency,  k,  for  which  there  is  experimental  data  is  less  than  k  ,  where  It1 

c4  c  4 

is  the  lowest  value  of  k  for  which  there  exists  a  solution  corresponding 

_1 

to  mode  4.  Therefore,  only  two  modes  will  be  of  interest  for  k  <  k 

c4 

The  modes  of  interest  are  modes  1  and  2  and  will  be  called  predominant 
modes . 

4.4  Summary 

The  results  of  the  solution  of  the  determinantal  equation  have  been 

presented  in  this  chapter.  It  was  seen  that  only  modes  1,  2  and  4  are 

1 

physically  admissable.  Mode  4  exists  only  for  k  >  k  .  Therefore,  for 
1  c4 
k  <  k  the  current  will  be  approximated  by  solutions  from  modes  1  and 


2. 


112 


5 .  SOURCE  PROBLEMS 


5.1  Introduction 

In  contrast  with  the  previous  chapters  that  considered  only  an  infinite 
tape  helix  without  a  source,  this  chapter  is  concerned  with  helices  with 
a  source.  The  source  will  be  of  small  extent  located  at  the  origin.  A 
description  for  the  source  is 


E 

11 


!<?' 


^o 


!4>i 


(5.1) 


where  $  is  much  less  than  one  radian, 
o 

taken  as 


In  the  limit  E  could  be 


En  =  6(<t>>  (5.2) 

Not  only  Will  the  infinite  helix  be  studied  but  also  a  finite  helix 

symmetrical  about  the  origin. 

5.2  An  Infinite  Helix  with  a  Source 

In  this  section  an  infinite  helix  with  a  source  as  described  in 

Equation  (5.1)  is  considered.  An  equation  for  the  amplitudes  of  the 

7 

free  modes  is  found  following  Sensiper  .  ,  In  a  range  of  frequencies 

of  interest  there  are  two  modes.  The  solution  for  the  relative  amplitudes 
is  shown  graphically. 


To  derive  the  expression  for  the  amplitudes  one  first  writes  the 
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If  the  integrand  is  set  equal  to  zero  an  expression  for  1(B)  is 
obtained.  Now  if  the  inverse  transform  of  1(B)  is  taken  the  result  is 


B4> 


I(4>)  =  I  e 
-00 


J 


.  E  sin  ,  2  1//2 

-i—  2  -169  .  ,  2  -  . 

_iB<{>  v  - b -  6  4ir  aA  (a  +  p  ) 


oo 


dB 


(5.6) 


I 


-iB <*>■*■  J  n2  2 r  2  ,,  ,1.  -2,  e"lkR  .,1 

e  <r-B  +k  [a  cos (9-9  )+p  j  — - —  d<$> 


R 


The  denominator  of  the  integrand  of  Equation  (5.6)  is  recognized  as  being 

proportional  to  E, ,  for  a  current  I  =  I  -iB<j>'*’.  The  electric  field  E, ,  for 
11  oe  11 

the  exponential  current  is  given  by  substituting  for  I  and  —  ,  in  Equation 

dq> 


(2.32)  viz. 


-i(jp 


'11  22  _2  1/2 

4irk  (a  +p  )  -oo 


00 

/ 


■>2,2.  2 


-B  -fk  [a  cos(4>-4>  )  +  p  ]>e 


_,1  ~ikR 

2 - d*1 


(5.7) 


Now  the  electric  field  may  also  be  found  from  Equation  (2.23) 


'll 


-i  B  — 

2p  -iB<^  .  2a, 

le _  I0ve _ sin^R 

4ir  k  £  p 


<2A  IB2  -  k 


4  [B 

o 


sin 


2* 


2  „  A  4  A  -2A 

k  ctn2^  _1  +1  ° 


2A„ 


(5.8) 


which  was  the  electric  field  found  in  the  derivation  of  the  determinantal 
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equation  by  matching  boundary  conditions. 

If  Equations  (5.7)  and  (5.8)  are  set  equal  to  each  other,  one  may 
substitute  for  the  denominator  of  the  integrand  of  Equation  (5.6)  the 
series  expression,  whose  zeros  have  been  found.  The  expression  for  I(4>) 
is 


T(<t» 


iE<t> 

o 


TT 


k  f-  p 

27 

tsin  + 


B<l> 

o 


"O'" 


6 

iB 

2p  -iBt 

e 


2 

k  ctn 


2a, 


dB 


(5.9) 


Let 


D(6) 


2 

k  c  t  n2'P 


(5 . 10) 


The  expression  D(B)  when  set  equal  to  zero  is  the  det erminantal 
equation.  The  roots  of  the  det erminantal  equation  are  poles  of  the 
integrand  in  Equation  (5.9).  The  integral  in  Equation  (5,9)  may  be 
evaluated  by  using  the  residue  theorem  to  find  the  current  in  terms  of 
the  free  modes  The  result  is 


4ttE4>  k  *  p 
o 


L 

1=1,2 


dt  D(B) 


B=B 

1 


K<!>)  = 


v  sin^  41 


(5,11) 
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where  use  was  made  of  the  fact, 


sxn  B.  — 


o 

”2 


1 


for  narrow  tapes  with  a  source  of  small  extent.  The  summation  in 
Equation  (5.11)  is  taken  over  the  two  predominant  modes. 

Let  the  current  on  the  tape  be 

-iB  4>  -iB  <j> 

I (4>)  =  IQ[e  +  a  e  ]  <j>  >  0  (5.12) 


A  calculation  for  the  magnitude  of  a,  |a|,  was  made  for  a  helix  with 
41  =  12.6°  and  §  =  .15,  The  results  are  shown  in  Figure  5.1.  Note  that 
as  the  frequency  k  is  increased  the  relative  amplitude  of  the  second 
mode  is  increased. 

The  complex  valued  derivatives  in  Equation  (5.11)  were  evaluated 
numerically  by  taking  first  an  increment  in  the  Br  direction  and  finding 


AD  AD 

r  i 

as  +1 

r  r 


(5.13) 


where  D  =  D  +  D. .  Second  an  increment  in  the  iB .  direction  was  taken 
r  1  1 


and  the  result 
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AD  AD 

_ i  r 

£6.  _1  m. 

1  l 

compared  to  Equation  (5.14)  to  determine  the  derivative.  This  is,  of 
course,  doing  nothing  more  than  performing  the  difference  quotient  along 
two  convenient  rays.  The  function  D(B)  is  analytic  in  a  region  containing 
the  zero  under  consideration,  and  the  derivative  is  unique.  The  only 
reason  that  the  difference  quotient  was  evaluated  along  two  different 
rays  was  to  have  a  numerical  check  on  the  calculation. 

5.3  A  Finite  Helix  with  Source 
5.3.1  Introduction 

This  section  will  be  concerned  with  relative  amplitudes  of  the 

two  predominant  current  modes  that  exist  on  a  finite  helix.  The  helix 

is  symmetrical  about  the  origin  -  <j>^<  ^  <  4>o  and  is  fed  at  the  origin. 

The  two  currents,  I  (  B  )  and  I  (B  ),  will  be  chosen  such  that  each 

current  will  be  zero  at  the  ends  of  the  helix,  <f>  =  +  <j*  .  This  boundary 

condition  excludes  consideration  of  mode  conversion  at  4>  =  +  4>  i,e., 

-  o  ’ 

no  energy  is  taken  from  mode  1  and  placed  in  mode  2  or  conversely 
by  the  end  discontinuity. 

The  currents  that  satisfy  the  boundary  condition  are 

Ix  =  sinh  Yx  (<j>o  -  I4>!)  (5.15) 


(5.14) 


and 


Ip  =  sin  B2(4>o~  I4>l) 


(5.16) 
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where  \  and  are  the  propagation  and  phase  constants  for  modes  one  and 
two  respectively,  which  were  found  by  solving  the  determinantal  equation 
for  an  infinite  helix  with  the  same  pitch  angle,  tape  width  and  radius  as 
the  finite  helix.  The  total  current  on  the  helix  is 

1  =  xi  +  AI  (5,17) 

where  A  is  a  complex  constant  which  relates  the  relative  amplitude  of  the 
two  modes. 

If  one  had  a  function  of  the  total  current  which  was  stationary  with 
respect  to  the  current  then  it  would  be  possible  to  find  the  complex 
constant  A  by  a  variational  technique.  This  is  how  A  is  calculated.  The 
relative  amplitude  of  the  two  modes  is  then  studied  as  a  function  of 
frequency  for  a  finite  helix  of  different  lengths. 

5,3.2  The  Variational  Expression 

A  variational  expression  that  is  stationary  with  respect  to  current 

30 

for  a  symmetric  helix  was  found  by  Tang 

The  integral  representation,  Equation  (2.31),  for  the  electric  field 
was  found  in  Chapter  2.  Equation  (2,31)  will  be  repeated  here  for 
convenience,  viz. 


E 


11 


-icjp. 


2,  2 

4trk  (a  + 


_2  1/2 
P  ) 


2  2  1  -2  e_:LkR  1 

+  k  [a  cos(.<(>-4>  )-*«p  ]  — - — 

R 


(5.18) 
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with 


2  cb-di1  n2 

R  =  [2a  sin  ?=|-] 


6  2 


+  [PC'M1)]2  »■ 


The  above  R  is  linearized  so  that  the  kernel  in  Equation  (5,18)  is  symmetric 
in  <)>  -  4*1,  i.e.  K($  -  =  K(<j>1  -  <j>) ,  with 


K  (<M  >  = 


-10/(1 


„  „  2  1/2 
A  ,2,  2  ~  \ 

417k  (a  -t-p  ) 


+  k2f a2cos(4-<t)1)+P  ]j 


-ikR 

e 

iT~ 


The  approximation  for  the  linearized  R  is  good  for  large  R  since  6 
is  compared  to  p ( 4*1 )  and  is  good  for  very  small  R  since  4  p ( 4^- 4*1 )  is 
compared  with  6. 

If  Equation  (5.18)  is  written  symbolically  as  E  =  K  I 
and  if 


<A,  B> 


ABd<j> 


where  the  helix  extends  from  - 
is  the  input  impedance 


<t> 

o 


+ 

o 


then  the  stationary  function 


Z. 

in 


<I,KI> 

i2(o) 


(5.19) 


in 

found  by  Tang  . 

The  stationary  function,  Equation  (5.19),  is  stationary  with 


respect  to  current,  i.e., 
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5  z 

- i-n  =  0 

6  I 

It  follows  from  Equations  (5.17)  and  (5.20)  that 


(5  20) 


6  Z 
_ in 

6  A 


=  0 


(5.21) 


If  the  indicated  operation  in  Equation  (5.21)  is  performed,  the  result 
for  the  complex  constant  A  is 


n  -  y°>  A2 

'  v0>  ST2 


(5,22) 


where 


I  (<j>) 
1 


( (t1) 


Kt^1)  d^cj)1 


(5  23) 


1^0)  -  sinh 


and 


I2(0)  =  310  Vo 

To  evaluate  the  complex  constant  A  it  is  necessary  to  evaluate  the  integrals 


in  Equation  (5.23), 
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5.3.3  Integration  of 


ij 


The  kernel  in  Equation  (5.23)  may  be  written  as 


K^1)  =  - 


„  ,  2 ,  2  -  .1/2  9<j>' 

4irk  (a  +p  ) 


d  d  9  i  t 

i  s»4> +  k  *<<m> V  oofr-*) 


(5.24) 


where 


f(u)  =  a  cos  u  4-  p 


-lkR 


G(u)  = 


with 


If 


2  r„  u,2  .2  ,6. 2 

R  =  [2a  sin  -]  +  (pu)  4-  (-) 


ij 


2  1/2 

4irk  (a  +-p  )  Q 
-iU|i  ij 


then 


l 


4>o 

1 — 

1  d<fn.(4>) 

*  (h 

i v+> 

8  d  9  1 

o 

o 

—  — * 

G(<j>-cj>1)d<j>1 

(5.25) 


The  inner  integral  is  separated  into  two  integrals,  viz. 
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0 

_ 

j  v*1* 

■t 

9  9 

8  4)1  ^ 

+  k2f(43-4>1) 

G(<j^4)1)  d4>1 


f 


+  f 

'0 


r 


9  9  ,  2,  ,  ,1 


9  cjs' 


1  ffej> 


+  k  (<M  ) 


g(4*“4^)  d4>1 


In  the  first  integral  above  the  change  of  variable  4*1  =  -  4*1  is  made  and 

in  addition  use  is  made  of  I  (-4*1)  =  I  (4*1)  to  yield 

J  j 


ij  =  J  d^I.(4>)  <  J  ij(4>1)  ~  [g^1)  -  g(<4^-4>1)  ]  d4>X 

-<f>  W>  9<^ 


(5.26) 


+  j  k2[f(<k4>1)  g(4h-4>1)  +  f(4>-4>1)  g(4>-4>1)]  d4>1'“ 


The  function  of  <j>  inside  the  braces, is  an  even  function  as  is  the 
current  I  (4>)  .  This  simplifies  Equation  (5.26)  to 


\  =  J  d4*io(4>)  <  J  I. OP1)  [g(<M>x)  -  g(4*-4>^) ]  V 

"o  (Vo  9t^ 

t,  1 

+  r  Ij c^1)  K2[f (cH^1)  g(4h-4>1)  +  f (4>-4>1)  G(<j^4>1)]  d4>xy  (5.27) 


The  first  integral  in  Equation  (5.27)  may  be  integrated  by  parts  twice, 
once  with  respect  to  4*  and  once  with  respect  to  4^  which  gives 


dl.  (4) 
1 


dl .(41) 
J 

d<4 


[G(4t41)  -  G(4‘-4>X)  ]  d^d^1 


(5,28) 


where  use  has  been  made  of 


1.(4)  =  1.(4)  =0  and 
l  o  jo 


G(41)  -  G  (-41)  =  0 


If  Equation  (5.28)  is  substituted  into  Equation  (5.27)  the  result  is 


[G(4r41)  -  g(4-41)]  d4d4X 


(5.29) 


+ 


ij(4x)  r(4)  k2[f (4f41)  g(4+-41)  +f (4-41)  G(<t^41).]d4d41 


It  will  be  reiterated  that  Equation  (5.29)  holds  only  when  the  following 
are  true: 


1)  I. .1.  are  even  functions  and  vanish  at  4 
ij  J  o 
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2)  R  is  linearized,  i.e.,  R(4)-<i>S  =  R(4>^_tl5)  • 

The  integrals  in  Equation  (5.29)  are  two  dimensional  integrals  and  need 
to  be  integrated  once  analytically  to  be  suitable  for  numerical  computation. 
The  reason  for  this  is  that  the  integrands  are  rapidly  varying  functions 
of  the  variable  of  integration  and  numerical  integration  requires  either 
a  very  small  step  size  or  a  variable  step  size  technique.  Either  of  the 
two  techniques  is  prohibitive  in  the  amount  of  machine  time  required,  for 
two  dimensional  integrals. 

However,  one  is  indeed  fortunate  that  after  a  change  of  variable  the 
integrals  may  be  integrated  analytically  once.  The  resulting  one 
dimensional  integrals  are  then  numerically  integrated  by  a  variable  step 
size  technique. 

Let 


u  =  (Jh^1 
v  _ 


(5. .30) 


Recall  that 


I  (40  =  sinh  y  (<(>  -4>) 

1  To 

i2(4>1)  =  sin  B2(4>o-4>1) 


dl1(4>) 

d<J> 

dI2(4>X) 

d4>x 


=  -  Y,  cosh  Y,  (4>  -40 

1  1  O 

=  -  B  cos  B2(45o_4)1) 


If  the  change  of  variable  in  Equation  (5.30)  is  made  then 


II (4>)  I1(4>1)  =  ~  [cosh  -  cosh  YjV] 


(5.31) 
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dl1(4>) 


d  I  x  (  ) 

d^1 


[  cosh  Y1(24>o-u)  +  cosh  -y^v] 


I  19)  l2l«P 


-  cosh  V„(24)  -u) 
'2  o 


cosh  Ygv  + 


sinh 


Y0(24>  -u) 

‘2  o 


sinh 


dl  (<|>)  d I  ( 4s1 ) 

S1- 


=  — —  [cosh  Y3(24>o-u)  cosh  Y2v  ~ 


sinh  Yg (240~u) 


sinh  Y„ 
’2 


+  cosh  Y  (2<f>  -u)  cosh  yv  -  sinh  y„(24  -u)  sinh  Y~vl 
2  o  O  2  0  6 


I2(4>)  Ig^1)  =  \  '(cos  B2v  -  cos  S2(24o-u)^ 


dl  (<j>)  dl  (41) 

A  A 


d4 


d<j) 


— -  [cos  B„  (24  -u)  +  cos  B  v] 
2  2  o  2  J 


where 


^2  = 


Y  B 

Yl+i  2 


(5.32) 


v 


(5.33) 


v 


(5.34) 


(5.35) 


(5.36) 


2 
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and 


iB2 


Since  the  assumed  currents  are  exponential  in  form  and,  of  course, 
have  the  property  e^X+^=eXe^,  the  products  of  the  currents  and  their 
derivatives  separate  into  products  of  functions  of  the  form  U(u)  V(v), 
where  U  is  a  function  of  u  alone  and  V  is  a  function  of  v  alone.  This 
particular  property  of  the  currents  permits  the  two  dimensional  integrals, 
Equation  (5.29)  to  be  integrated  once  analytically. 

The  integrals,  Equation  (5.29),  are  now  of  the  form 


(<H4>^)  F^cjj-qT)  d<j>d4>1  (5.37) 


The  Jacobian  of  the  transformation  given  in  Equation  (5  30)  is  —  .  The 
integral  in  Equation  (5.37)  may  then  be  written  as 


5  k  ■  //  F1(u)  Fg(v)  dudv  (5.38) 

R 

uv 


where  R  is  shown  in  Figure  5.2. 
uv 
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Figure  5.2 

If  the  integration  with  respect  to  u  is  to  be  integrated  first,  i,e., 
able  to  be  integrated  analytically,  Equation  (5  38)  is  written  as 


1 

2 


F2(v)  +  F2 (-v) 


dv 


i 


2<f>  -V 
o 


F^ (u) du 


(5.39) 


On  the  other  hand  if  the  integration  with  respect  to  v  is  to  be  performed 
first,  Equation  (5.38)  is  written 


1. 

2 


u 


f2(v) 


dv 


(u)  du 


L 


2<f>  -u 
-  o 


Fg(v) dv 


-2<t>  +u 
o 


(5  o  40) 
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If  the  currents  and  their  derivatives,  Equations  (5,31)  through  (5,36), 
are  substituted  into  Equations  (2.39)  and  (3, .40),  and  the  integrations 
first  are  performed,  the  results  for  the  integrals,  Equation  (5.29), 
are  obtained  after  considerable  analytical  integration.  The  results  are 


11 


S’  f 

/°  2  2 

2[V1  *■  k  f(u)]  [u  cosh  ^(24^-u)  -2  (<j>  -u)  cosh  Y-^u] 


+[ V12-k2f (u) ] 


2  sinh  y^u  -  sinh  Y^^^-u) 


G ( u)  du 


2<t> 

r 

Js 


2  ,  2 


+  I  /[^  +k  f  (u)  )  (2^-u)  cosh  Yl(2tj>o-u) 


9  o  sinh  Yi  (2<f>  -u) 

4-  -k  f  ( u)  ]  - — - - - f  G (u)  du 


(5.41) 
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4> 
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G(u) 


bi^g-i  k2f(u)] 


sinh  Y3U 

- — - cosh  [  Y„  ( 2  41  -u) 

Y  *2  To 


-  cosh  Y„u 
'2 


sinh  [  Yg  (2^-u)  ]  -  sinh  Y3u 


2  JSltlh  V1 
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sinh  [ v  (2<f>  -u)]  -  sinh  v  u 
2  o  ' 


cosh  "Yg11 


+  |  G(u)  m^-ik  f(u)J 


sinh  y, -u) 

- — - cosh  [  V  ('24’  -u)  ] 

v  2  o 

v3 


(5.42) 


sinh  Y-(24)  -u)  | 

+  [Y,  B  +  i  k2f(u)]  - - - - -  cosh  [V3<2<i)0-u)U  du 

12  Y2 


Zo  2  sin  B  u  -  sin  B  (2$  -u) 

22  f  [B/  +  A<»» - S— — - 


+  [B  2-k2f(u)]  fu  cos  B  (24>o-u)  -2  (<j>Q-u)  cos  B2u]  K  G(u)  du 


r9of  o  o  sin  B _(2«t»  -u) 

+  f  J[S22  +  k2  f(u>] - §——  +  [B2  -k 


+  [B  -k  f(u)]  (24>o-u)  cos  S 2 ( 2 <t>o — u.)r  G(u) du 


(5.43) 


f(u)  =  a  cos  u  +  p 


e"lkR  2  ,  .  u t2  . 2  A 

G(u)  =  — - —  R  =  [2a  sin  -j  +  (pu)  f  (-) 


V  B 
Yl+i  2 

Y2  2  J 


v  B 

'1-i  i 


Y3  ”  2" 


where 
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The  integrals  are  now  in  a  form  suitable  for  numerical  evaluation  and 
hence  the  complex  constant  A  may  be  found  from  Equation  (.5,22)  which 
is  repeated  here  for  convenience  after  the  numerator  and  denominator  on 
the  right  hand  side  are  multiplied  by  a  constant. 


A 


ywi^i  -  yo>^2 

-  V°>*£ 


(5  44) 


with 


and 


I,  (0)  =  sinh  v  <j> 

1  1  o 


I2(0)  =  sin  B2<t>o 

5.3.4  Results  of  the  Numerical  Integration 

The  integrals  in  Equation  (5.44)  were  evaluated  using  a  variable  step 
size  technique  and  the  value  of  A  was  determined  for  helices  with  a  total 
of  8,  12,  and  16  turns ,  All  helices  have  a  pitch  angle  41  =  12  6  and  a 
tape  width  ®  =  .15. 

One  is  interested  in  how  much  of  the  second  mode  is  launched  as 

compared  to  the  first  mode.  For  this  reason  the  results  for  A  are  not 

given,  but  the  values  for  |a|  =  are  given.  The  magnitude 

S 1 n  Xl  o 

of  a,  | a  |,  approximately  equals  the  ratio  of  the  amplitude  of  the  first 

~2Vo 

mode.  The  approximation  arises  from  ignoring  the  term  e  as 

compared  to  unity.  The  approximation  is  better  than  2 jo  on  the  worst  case 
and  often  better  than  .05 'jo. 
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In  Figure  5.3  I  a |  is  plotted  as  a  function  of  frequency  for  the 
three  different  length  helices.  It  is  to  be  noted  as  frequency  is 
increased  the  general  trend  is  for  an  increase  in  second  mode.  This 
functional  variation  is  identical  to  what  was  found  for  the  infinite 
helix  as  shown  in  Figure  2.1. 

The  oscillations  in  the  amplitude  of  a  may  be  caused  by  the  variation 

of  input  impedance  of  the  helix.  The  minimums  of  j  a  J  occur  near  values 

of  k  which  correspond  to  values  of  B  „  that  make  the  helix  resonant, 

r2 

The  magnitude  of  the  total  current  on  a  helix  of  12  turns  for  two 
different  frequencies  is  shown  in  Figures  5.4  and  5.5.  At  the  lower 
frequency,  k  =  .17,  the  attenuation  for  the  first  mode  was  smaller  than 
at  the  higher  frequency,  k  =  .29.  Consequently,  the  total  current  does 
not  decay  as  fast  at  the  lower  frequency.  In  addition,  the  amplitude  of 
the  second  mode  is  less  at  the  lower  frequency.'  Therefore,  the  second 
mode  does  not  give  the  pronounced  standing  wave  as  at  the  higher  frequency. 
Figures  5.6a,  6b,  and  6c  show  that  k  -  B  diagram  in  the  region  of  interest 
in  this  chapter. 

5 . 4  Summary 

Two  source  problems  were  investigated  in  this  chapter,  one  an  infinite 
helix  and  the  other  a  finite  helix.  Both  problems  were  similar  in  that  the 
helices  were  fed  at  the  origin  with  sources  of  small  extent,  essentially  6  - 
sources.  In  both  problems  the  object  was  to  find  the  relationship  between 
the  amplitudes  of  the  two  predominant  modes  constituting  the  current  on  the 
helix.  The  phase  constants  determined  by  solving  the  determinantal  equation 
in  previous  chapters  were  used  in  finding  the  relative  amplitudes  of  the 


two  modes . 


1.0 
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k  — ► 

Figure  5.3.  Relative  Amplitude  of  the  Two  Modes  on  a  Finite 
'  Helix  as  a  Function  of  Frequency. 
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Figure  5.6a.  The  k  -  8  diagram  for  the  tape  helix  near 

k  , ;  Moderi;  ^  =  12.6°,  a  =  .335. 
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Figure  5.6b.  The  k  —  B.  diagram  for  the  tape  helix  near 
kcl;  Mode1!;  +  =  12.6°,  a  -  .335. 
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Figure  5.6c. 


k  -* 

The  k  -  B  diagram  for  the  tape  helix  near 
k  ;  Moder2;  +  =  12.6°,  a  =  .335. 

C  X. 
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The  relative  amplitude  of  the  two  modes  in  the  case  of  the  infinite 
helix  was  found  by  evaluating  the  contribution  made  by  the  leaky  wave  poles 
to  the  inverse  Fourier  transform  of  the  source  problem.  The  leaky  wave 
poles  are  the  solutions  to  the  determinantal  equation. 

The  relative  amplitude  of  the  two  modes  in  the  case  of  the  finite 
helix  was  found  by  using  a  variational  formulation  to  find  the  relative 
amplitude  to  give  the  best  input  impedance  in  terms  of  the  two  current 
waves  which  correspond  to  solutions  of  the  determinantal  equation.. 

It  was  found  that  the  second  mode  increased  in  amplitude  as  the 
frequency  was  increased  for  both  problems.  In  addition,  when  the  helix 
was  truncated,  there  was  an  end  effect  and  the  oscillations  of  the 
relative  amplitude  agree  well  with  the  expected  variations  of  the  input 


impedance  for  the  second  mode. 
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6.  COMPARISON  WITH  EXPERIMENTS 


6.1  Marsh's  Experiment 

6.1.1  Introduction 

In  1950  Marsh*-^  reported  that  he  had  measured  the  current  distri¬ 
bution  of  a  six  turn  helix  fed  against  a  ground  plane.  In  addition  to 
measuring  the  current  he  was  able  to  empirically  fit  the  results  by 
assuming  two  current  modes  each  of  which  was  reflected  from  the  end. 

For  the  empirical  results  he  chose  the  current  as 

-(a  +iS  )z  -iB  z  (a  +iB  )z  iB  z 
_  JL  i.  £  ,  JL  -L  £ 

I  =  e  +  ae  +  be  +  ce 

The  real  constants  ,  6^,  and  S2  and  the  complex  constants  a,  b, 
and  c  were  chosen  by  Marsh  to  fit  his  measured  data.  He  obtained  excellent 
agreement . 

6.1.2  Phase  Constants 

The  phase  constants  calculated  by  solving  the  determinantal  equation 
agree  well  with  Marsh’s  empirically  determined  phase  constants  as  shown 
in  TABLE  I. 

TABLE  I 


“l 

S1 

B2 

k 

Marsh 

Calculated 

Marsh 

Calculated 

Marsh 

Calculated 

.1488 

.0078 

0 

.681 

.69 

1.193 

1.15 

.1734 

.0636 

.069 

.795 

.808 

1.173 

1.175 

.2115 

.1548 

.1065 

.969 

.875 

1.356 

1.212 

141 


The  calculated  phase  constants  in  TABLE  I  correspond  to  tape  helix 
parameters,  41  =  12.6°,  a  =  4,305  cm  and  6  =  .15a.  The  empirically  deter¬ 
mined  phase  constants  by  Marsh  correspond  to  parameters  41  =  12.6°  and 
a  =  4.305  cm.  The  wire  diameter  is  not  reported  by  Marsh  and  a  tape 
width  was  chosen  to  be  approximately  equal  to  the  wire  diameter  inferred 
from  a  photograph  in  Marsh's  report.  For  Marsh  used  Sq  corresponding 
to  a  wave  traveling  down  the  wire  with  the  velocity  of  light. 

6.1.3  Amplitude  Constants 

The  solution  for  the  helix  fed  against  the  ground  plane  was  not 
attempted.  However,  it  is  interesting  to  compare  the  magnitude  of  the 
amplitude  coefficients  obtained  for:  1)  The  infinite  helix  fed  at  origin; 
2)  the  finite  helix  fed  at  the  origin,  and  3)  Marsh's  empirical  results 
for  the  helix  fed  against  a  ground  plane.  The  comparison  is  shown  in 
Figure  6.1.  Note  that,  as  k  is  increased,  one  in  general  expects  the 
second  mode  amplitude  to  increase  relative  to  the  first  mode  amplitude. 
Marsh,  in  empirically  determining  the  amplitude  coefficients,  did  not 
insist  on  zero  current  for  each  mode  or  total  at  the  end  of  the  helix. 

He  permitted  energy  propagation  down  the  helix  in  mode  2  to  be 
converted  into  mode  1  to  be  reflected.  This  is  reasonable  since  at  the 
discontinuity  the  waves  are  'Relaunched"  and  mode  1  is  easier  to  excite 
than  mode  2  as  shown  by  the  calculations  of  Chapter  5.  Since  only  a 
portion  of  the  mode  2  current  returns  this  may  permit  a  larger  mode  2 
current  to  be  launched  and  hence  a  larger  relative  amplitude,  a.  In 
Chapter  5,  where  source  problems  were  investigated,  each  mode  current  was 
forced  to  be  zero  at  the  end  of  the  helix. 
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k  — ► 

Figure  6.1.  Relative  amplitude  of  the  two  modes  compared  for 
a)  infinite  helix,  b)  finite  helix  and  c.)  Marsh's 
empirical  work. 
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6.1.4  Summary 

Marsh  was  able  to  empirically  fit  the  measured  current  distribution 
on  a  six  turn  helix  fed  against  a  ground  plane  by  assuming  two  current 
waves  each  reflected  at  the  end.  The  phase  constants  calculated  in  this 
report  agree  well  with  Marsh's  empirically  determined  phase  constants. 

In  addition  the  amplitude  coefficients  calculated  m  Chapter  5  behave 
similarly  to  Marsh's  coefficients. 

6.2  Other  Experiments 

31  12 

Both  McClelland;,  on  the  conical  equiangular  spiral,  and  Patton, 

on  the  bifilar  helix  fed  out  of  phase,  have  performed  experiments  which 

support  the  use  of  the  approximation  for  the  current  by  only  the  two 

modes.  The  bifilar  helix  fed  out  of  phase  has  a  determinantal  equation 

which  is  related  to  the  determinantal  equation  for  the  monofilar  helix 

studied  in  this  report.  The  determinantal  equation  for  the  bifilar 

helix  has  only  the  terms  corresponding  to  n  odd  as  coefficients  to  D 

as  compared  to  the  determinantal  equation  for  the  monofilar  helix.  One 

expects  similar  behavior  of  the  solution  near  =  2  but  certainly  different 

solutions  near  B  =2  since  the  terms  containing  t  are  missing.  The 
r  b  -2 

conical  equiangular  spiral  is  the  log-periodic  version  of  the  bifilar 
helix.  Both  works,  McClelland's  and  Patton's,  support  the  calculations 
reported  here  in  the  following  manner: 

1)  As  the  frequency  is  increased  above  a  value  corresponding  to 
k f  the  current  decays,  indicating  radiation.  The  higher  k 

the  more  attenuation,  indicating  B.^  increasing  as  k  is  increased. 

2)  As  k  is  increased  still  further  the  current  has  a  small 


standing  wave  component  occurring  indicating  that  mode  2  is  present. 
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The  foregoing  experimental  results  are  exactly  what  would  be  predicted 
from  the  solution  of  the  detenn inantal  equation  and  the  corresponding  study 
of  the  source  problem. 


I 
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7.  SUMMARY 


7 , 1  Original  Work  Done 

Sensiper  formulated  and  solved  the  determinantal  equation  for  the 
tape  helix  treating  the  phase  constant,  B,  as  a  real  variable  but 
Sensiper  and  Pierce  doubted  the  existence  of  complex-valued  solutions 
for  the  determinantal  equation.  This  present  work  reports  the  solution 
of  the  determinantal  equation  for  various  values  of  the  parameters  of 
the  tape  helix,  where  the  phase  constant  is  complex-valued.  The 
detei minantal  equation  was  interpreted  as  a  complex-valued  equation  where 
the  arguments  of  the  IK  products  are  shown  to  be  continuous  functions 
of  k  and  B.  In  determining  how  the  various  arguments  of  the  IK  products 
should  be  interpreted  it  was  found  useful  to  study  two  simplified  equations 
which  exhibited  solutions,  in  certain  regions  of  B^,  similar  in  character 
to  the  solutions  for  the  determinantal  equal  ion. 

Two  source  problems  were  investigated,  the  first  an  infinite  helix, 
the  second  a  finite  helix,  both  fed  at  the  origin.  The  amplitudes  of 
the  free  modes  representing  the  current  on  the  infinite  helix  were  found 

7 

following  Sensiper  but  with  the  phase  constant,  B,  as  a  complex 
variable.  On  the  problem  of  a  finite  helix  a  variational  technique  was 
used  to  find  the  relative  amplitude  of  the  two  modes  representing  the 
current  which  yields  the  best  input  impedance  for  a  current  approximated 

by  two  modes  which  correspond  to  solutions  of  the  determinantal  equation. 

30 

Tang  used  the  variational  technique  to  find  the  input  impedance  of 
the  helical  antenna  (finite  helix)  assuming  the  current  to  be  the  sum  of 
two  components,  the  first,  was  sinusoidal  and  corresponded  to  waves 
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traveling  at  the  velocity  of  light;  the  second  was  also  sinusoidal  but 
with  the  phase  constant  twice  the  value  corresponding  to  phase  constant 
for  a  wave  traveling  at  the  velocity  of  light. 

The  variational  technique  was  used  in  the  present  report  to  find 
the  relative  amplitude  of  the  two  assumed  modes  which  were  used  to 
represent  the  current  on  the  finite  helix.  The  phase  constants 
corresponding  to  the  assumed  modes  are  the  complex-valued  solutions 
to  the  determinantal  equation. 

7.2  Conclusions 

This  report  indicates  that  it  is  possible  to  approximate  the 
current  on  a  helix  with  current  waves  whose  complex-valued  phase  constants 
are  found  solving  the  determinantal  equation.  The  value  and  the  behavior 
of  the  roots  as  a  function  of  the  helix  parameters  of  the  determinantal 
equation  is  Just  that  behavior  which  causes  the  corresponding  current 
to  agree  with  experimental  results. 

“308 

The  solution  of  the  determinantal  equation  for  a  =  10  an 

extremely  narrow  tape,  corresponds  to  the  asymptotic  solution  obtained 

by  Mittra^1  for  the  limiting  case  where  the  tape  width  tends  to  zero. 

As  the  tape  is  made  wider  but  still  very  narrow,  e.g.,  corresponding  to 
-30 

a  =  10  ,  the  phase  constant  becomes  complex  valued  and  in  addition 

the  real  part  of  the  phase  constant  deviates  from  the  line  k  =  6  sin  4* 

r 

representing  a  wave  traveling  down  the  tape  with  the  velocity  of  light. 

The  wider  the  tape  the  greater  is  the  value  of  8^,  the  imaginary  part 
of  the  phase  constant.  Also  the  wider  the  tape  the  greater  is  the 
deviation  of  the  real  part  of  the  phase  constant  from  the  line 


corresponding  to  the  velocity  of  light. 
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For  the  mode  1  solution  the  imaginary  part  of  the  phase  constant, 
increases  rapidly  as  the  frequency  is  increased  above  k  =  k^. 

As  the  pitch  angle,  4*,  is  increased  the  real  part  of  the  phase 
constant  tends  toward  the  line  corresponding  to  the  velocity  of  light. 

Also  as  the  pitch  angle  is  inci.  eased  the  maximum  value  of  the  imaginary 
part  of  the  phase  constant,  decreases. 

The  amplitude  coefficients  calculated  in  connection  with  the  source 
problems  permit  the  calculation  of  the  total  current  on  the  helix. 
Corresponding  to  a  frequency  just  above  k  -  k  the  current  is 
predominantly  mode  1  type  current  and  shows  exponential  decay  away 
from  the  source.  For  a  higher  value  of  k  the  current  shows  the  relative 
increase  of  mode  2  type  which  is  typified  by  a  definite  standing  wave 
of  current  near  the  reflected  end  for  the  finite  helix, 

7,3  Further  work 

7.3.1  The  Bifilar  Helix 

The  determi nantal  equation  has  been  derived  for  the  bifilar  helix  when 
the  currents  are  assumed  balanced.  This  det erminantal  equation  has  only 
the  terms  on  the  right  hand  side  which  have  form  factors,  i.e.,  the  even 
terms  are  zero.  Although  simplified  Equation  (1),  Section  3.2,  gives  an 
insight  into  the  behavior  of  the  solution,  a  more  exact  solution  might  be 
found.  After  one  obtains  the  roots  for  the  infinite  bifilar  helix,  the 

source  problem  could  be  investigated  with  the  end  results  the  calculation 

1 2 

of  patterns.  These  could  be  checked  against  Patton's  calculated  and 
measured  results. 

7.3.2  _ A  Conducting  Cylinder  Inside  the  Helix 

Following  the  method  similar  to  that  used  in  this  report^  the 
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determinantal  equation  for  a  concentric  conducting  cylinder  inside 
the  helix  may  be  obtained.  The  behavior  of  the  solutions  could  then 
be  used  to  verify  the  preliminary  measurements  12  indicating  that 
the  radius  of  cylinder  has  small  effect  on  the  fields  outside  the 


helix. 
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Bolling  Air  Force  Base 
Washington  25,  D„  C. 

Office  of  Chief  Signal  Officer 
Engineering  &.  Technical  Division 
Attn:  SIGNET- 5 
Washington  25,  D.  C. 

Commanding  Officer 

U.  S.  Army  Electronics  R&D  Activity 
Attn:  SELWS-ED 

White  Sands,  Missile  Range,  New  Mexico 
Director 

Surveillance  Department 
Evans  Area 

Attn:  Technical  Document  Center 
Belman,  New  Jersey 

Commander 

U.  S,  Naval  Air  Test  Center 
Attn:  WST-54,  Antenna  Section 
Patuxent  River,  Maryland 

Material  Laboratory,  Code  932 
New  York  Naval  Shipyard 
Brooklyn  1,  New  York 

Director 

Naval  Research  Laboratory 
Attn:  Code  5200 
Washington  25,  D.  C. 


Director 

Air  University  Library 
Attn:  3T-AUL-59-30 
Maxwell  Air  Force  Base 
Alabama 

Commanding  General 

White  Sands  Missile  Range 

Attn:  ORDBS-OM-Tech  Library  RR-312 

New -Mexico. 

Commanding  Officer 

Diamond  Ordnance  Fuse  Laboratories 

Attn;  240 

Washington  25,  D.  C. 

Director 

U.  S.  Navy  Electronics  Laboratory 

Attn:  Library 

San  Diego  52,  California 

Adamt-Russell  Company 

200  Sixth  Street 

Attn:  Library  (Antenna  Section) 

Cambridge,  Massachusetts 

Aero  Geo  Astro 
Attn:  Security  Officer 
1200  Duke  Street 
Alexandria,  Virginia 

NASA  Goddard  Space  Flight  Center 
Attn:  Antenna  Section,  Code  523 
Greenbelt,  Maryland 

Airborne  Instruments  Labs.,  Inc. 
Attn:  Librarian  (Antenna  Section) 
Walt  Whitman  Road 
Melville,  L.I.,  New  York 

American  Electronic  Labs 
Box  552  (Antenna  Section) 

Lansdale,  Pennsylvania 

Andrew  Alfred  Consulting  Engineers 
Attn:  Librarian  (Antenna  Section) 
299  Atlantic  Avenue 
Boston  10,  Massachusetts 
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Amphenol-Borg  Electronic  Corporation 
Attn:  Librarian  (Antenna  Section) 
2801  Sa  25th  Avenue 
Broadview,  Illinois 

Bell  Aircraft  Corporation 
Attn:  Technical  Library 

(Antenna  Section) 

Buffalo  5,  New  York 

Bendix  Radio  Division  of 

Bendix  Aviation  Corporation 
Attn:  Technical  Library 
(For  Dept  .  462-4) 

Baltimore  4,  Maryland 

Antenna  Systems,  Inc. 

Manager  Documentation 
Grenier  Field 
Manchester,  New  Hampshire 

Boeing  Airplane  Company 
Aero  Space  Division 
Attn:  Technical  Library 

M/F  Antenna  &  Radomes  Unit 
Seattle,  Washington 

Boeing  Airplane  Company 
Attn:  Technical  Library 

M/F  Antenna  Systems  Staff  Unit 
Wichita,  Kansas 

Chance  Vought  Aircraft,  Inc. 

THRU:  Bu  AER  Representative 
Attn.  Technical  Library 

M/F  Antenna  Section 
P.  O.  Box  5907 
Dallas  22,  Texas 

Collins  Radio  Company 
Research  Division 
Attn:  Technical  Library 
Cedar  Rapids,  Iowa 

Convair 

Ft.  Worth  Division 
Attn;  Technical  Library  (Antenna 
Section) 

Grants  Lane 
Fort  Worth,  Texas 


Convair 

Attn .  Technical  Library 
(Antenna  Section) 

P.  0.  Box  1050 

San  Diego  12,  California 

Dalmo  Victor  Company 
Attn:  Technical  Library 
(Antenna  Section) 

1515  Industrial  Way 
Belmont,  California 

Dome  &  Margolin,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

30  Sylvester  Street 
Westbury,  L.  I.,  New  York 

Dynatronics,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

Orlando,  Florida 

Electronic  Communications,  Inc. 

Research  Division 
Attn:  Technical  Library 
1830  York  Road 
Timonium,  Maryland 

Fairchild  Engine  &  Airplane  Corporation 
Fairchild  Aircraft  &  Missiles  Division 
Attn  Technical  Library 
(Antenna  Section) 

Hagerstown  10,  Maryland 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
Attn-  Technical  Library 

M/F  Electronics  Division 
Atlanta  13,  Georgia 

General  Electric  Company 
Electronics  Laboratory 
Attn:  Technical  Library 
Electronics  Park 
Syracuse,  New  York 

General  Electronic  Labs.,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

18  Ames  Street 

Cambridge  42,  Massachusetts 
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General  Precision  Lab.,  Division  of 
General  Precision,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

63  Bedford  Road 
Pleasantville,  New  York 

Goodyear  Aircraft  Corporation 
Attn:  Technical  Library 
M/F  Dept  .  474 
1210  Massilon  Road 
Akron  15,  Ohio 

Granger  Associates 
Attn:  Technical  Library 
(Antenna  Section) 

974  Commercial  Street 
Palo  Alto,  California 

I-T-E  Circuit  Breaker  Company 
Special  Products  Division 
Attn:  Research  Library 
601  E,  Erie  Avenue 
Philadelphia  34,  Pennsylvania 

General  Mills  Electronics  Division 
Attn:  Dr,  H,  P.  Raabe 
2003  East  Hennepin  Avenue 
Minneapolis  13,  Minnesota 

Grumman  Aircraft  Engineering  Corporation 
Attn;  Technical  Library 

M/F  Avionics  Engineering 
Bethpage,  New  York 

The  Hallicraf ters  Company 
Attn:  Technical  Library 
(Antenna  Section) 

4401  W.  Fifth  Avenue 
Chicago  24,  Illinois 

Hoffman  Laboratories,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

Los  Angeles  7,  California 

John  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Springs,  Maryland 


Hughes  Aircraft  Corporation 
Attn:  Technical  Library 
(Antenna  Section) 

Florence  &  Teal  Street 
Culver  City,  California 

ITT  Laboratories 
Attn:  Technical  Library 
(Antenna  Section) 

500  Washington  Avenue 
Nutley  10,  New  Jersey 

U.  S.  Naval  Ordnance  Laboratory 
Attn:  Technical  Library 

Corona,  California 

Lincoln  Laboratories 

Massachusetts  Institute  of  Technology 
Attn:  Document  Room 
P .  O .  Box  7  3 

Lexington  73,  Massachusetts 

Litton  Industries 
Attn:  Technical  Library 
(Antenna  Section) 

4900  Calvert  Road 
College.  Park,  Maryland 

John  Hopkins  University 
Radiation  Laboratory 
Attn:  Library 
1315  St.  Paul  Street 
Baltimore  2,  Maryland 

Lockheed  Missile  &  Space  Division 
Attn;  Technical  Library  (M/F  Dept- 
58-40,  Plant  1,  Bldg.  130) 
Sunnyvale,  California 

The  Martin  Company 
Attn:  Technical  Library 

(Antenna  Section)  Mail  No.  T-38 
P.  0.  Box  179 
Denver  1,  Colorado 

The  Martin  Company 
Attn:  Technical  Library  (Antenna 
Section)  .  ,> 

Baltimore  3,  Maryland 
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The  Martin  Company 
Attn:  Technical  Library  (M/F 
Microwave  Laboratory) 

Box  5837 

Orlando,  Florida 

W.  L.  Maxson  Corporation 
Attn:  Technical  Library 
(Antenna  Section) 

460  West  34th  Street 
New  York  1,  New  York 

McDonnell  Aircraft  Corporation 
Attn:  Technical  Library 
(Antenna  Section) 

Box  516 

St.  Louis  66,  Missouri 
Melpar,  Inc. 

Attn:  Technical  Library 
(Antenna  Section) 

3000  Arlington  Blvd . 

Falls  Church,  Virginia 

University  of  Michigan 
Radiation  Laboratory 
Willow  Run 
201  Catherine  Street 
Ann  Arbor,  Michigan 

Mitre  Corporation 

Attn:  Technical  Library  (M/F  Electronic 
Warefare  Dept.  D-21) 

Middlesex  Turnpike 
Bedford,  Massachusetts 

Northeastern  University 
Attn:  Dodge  Library 
Boston  15,  Massachusetts 

North  American  Aviation,  Inc. 

Attn:  Technical  Library  (M/F 
Engineering  Dept.) 

4300  E.  Fifth  Avenue 
Columbus  16,  Ohio 

North  American  Aviation,  Inc. 

Attn:  Technical  Library 
(M/F  Dept.  56) 

International  Airport 
Los  Angeles,  California 


Northrop  Corporation 
NORAIR  Division 
1001  East  Broadway 

Attn:  Technical  Information  (3924-3) 
Hawthorne,  California 

Ohio  State  University  Research 
Foundation 

Attn:  Technical  Library 

(M/F  Antenna  Laboratory) 

1314  Kinnear  Road 
Columbus  12,  Ohio 

Ohio  State  University  Research 
Foundation 

Attn:  Dr.  C  „  H.  Walter 
1314  Kinnear  Road 
Columbus  12,  Ohio 

Philco  Corporation 

Government  &  Industrial  Division 

Attn:  Technical  Library 

(M/F  Antenna  Section) 

4700  Wissachickon  Avenue 
Philadelphia  44,  Pennsylvania 

Westinghouse  Electric  Corporation 

Air  Arms  Division 

Attn:  Librarian  (Antenna  Lab) 

P.  0.  Box  746 
Baltimore  3,  Maryland 

Wheeler  Laboratories 

Attn:  Librarian  (Antenna  Lab) 

Box  56 1 

Smithtown,  New  York 

Electrical  Engineering  Research 
Laboratory 
University  of  Texas 
Box  8026,  University  Station 
Austin,  Texas 

University  of  Michigan  Research 
Institute 

Electronic  Defense  Group 
Attn:  Dr.  J.  A.  M.  Lyons 
Ann  Arbor,  Michigan 

Radio  Corporation  of  America 
RCA  Laboratories  Division 
Attn:  Technical  Library 

(M/F  Antenna  Section) 
Princeton,  New  Jersey 
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Radiation,  Inc. 

Attn:  Technical  Library  (M/F) 

Antenna  Section 
Drawer  37 

Melbourne,  Florida 

Ramo-Wooldridge  Corporation 
Attn:  Librarian  (Antenna  Lab) 

Canoga  Park,  California 

Rand  Corporation 

Attn:  Librarian  (Antenna  Lab) 

1700  Main  Street 
Santa  Monica,  California 

Rantec  Corporation 

Attn:  Librarian  (Antenna  Lab) 

23999  Ventura  Blvd . 

Calabasas,  California 

Raytheon  Corporation 
Equipment  Division 
Library  -  J.  Portsch 
P.  O.  Box  520 
Waltham  54,  Massachusetts 

Republic  Aviation  Corporation 
Applied  Research  &  Development 
Division 

Attn:  Librarian  (Antenna  Lab) 
Farmingdale,  New  York 

Sanders  Associates 

Attn:  Librarian  (Antenna  Lab) 

95  Canal  Street 
Nashua,  New  Hampshire 

Southwest  Research  Institute 
Attn:  Librarian  (Antenna  Lab) 

8500  Culebra  Road 
San  Antonio,  Texas 

H.  R.  B.  Singer  Corporation 
Attn:  Librarian  (Antenna  Lab) 

State  College,  Pennsylvania 

Sperry  Microwave  Electronics  Company 
Attn:  Librarian  (Antenna  Lab) 

P.  0.  Box  1828 
Clearwater,  Florida 


Sperry  Gyroscope  Company 
Attn:  Librarian  (Antenna  Lab) 
Great  Neck,  L.  I.,  New  York 

Stanford  Electronic  Laboratory 
Attn:  Librarian  (Antenna  Lab) 
Stanford,  California 

Stanford  Research  Institute 
Attn:  Librarian  (Antenna  Lab) 
Menlo  Park,  California 

Sylvania  Electronic  System 
Attn:  Librarian  (M/F  Antenna  & 
Microwave  Lab) 

100  First  Street 
Waltham  54,  Massachusetts 

Sylvania  Electronic  System 
Attn:  Librarian  (Antenna  Lab) 

P.  0.  Box  188 

Mountain  View,  California 

Technical  Research  Group 
Attn:  Librarian  (Antenna  Section) 
2  Aerial  Way 
Syosset,  New  York 

Ling  Temco  Aircraft  Corporation 
Temco  Aircraft  Division 
Attn:  Librarian  (Antenna  Lab) 
Garland,  Texas 

Texas  Instruments,  Inc. 

Attn:  Librarian  (Antenna  Lab) 

6000  Lemmon  Avenue 
Dallas  9,  Texas 

A.  S.  Thomas,  Inc. 

Attn:  Librarian  (Antenna  Lab) 

355  Providence  Highway 
Westwood,  Massachusetts 

New  Mexico  State  University 
Head  Antenna  Department 
Physical  Science  Laboratory 
University  Park,  New  Mexico 
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Bell  Telephone  Laboratories,  Inc. 
Whippany,  New  Jersey 
Attn-  Technical  Reports  Librarian 
Room  2A-1.65 

Robert  G.  Hansen 
Aerospace  Corporation 
Box  95085 

Los  Angeles  45,  California 

Dr.  Richard  C,  Becker 
10829  Berkshire 
Westchester,  Illinois 

Dr.  W.  M.  Hall 
Raytheon  Company 
Surface  Padj.r  and  Navigation 
Oper at  ions 
Boston  Post  Road 
Wayland,  Massachusetts 

Dr.  Robert  L.  Carrel 
Collins  Radio  Corporation 
Antenna  Section 
Dallas,  Texas 

Dr.  A.  K.  Chat  ter jee 
Vice  Principal  and  Head  of  the  De¬ 
partment  of  Research 
Birla  Institute  of  Technology 
P.  0.  Mesra 

District-Ranchi  (Bihar)  India 

University  of  Dayton 

Research  Institute 

Attn-  Professor  Douglas  Hanneman 

300  College  Park 

Dayton,  Ohio 

Technische  Hochschule 
Attn:  H.  H.  Meinke 
Munich,  Germany 

NASA  Goddard  Space  Flight  Center 
Attn:  Antenna  Branch 
Mr.  Lantz 

Greenbelt,  Maryland 

Aeronautic  Division 

Ford  Motor  Company 

Ford  Road  -  Attn:  Mr.  J.  M.  Black 

Newport  Beach,  California 


Professor  A.  A.  Oliver 
Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,.  New  York 

U.  S.  Naval  Ordnance  Laboratory 
Attn:  Technical  Library 
Corona,  California 

Avco  Corporation 

Research  and  Advanced  Development 
Division 

Attn  Research  Library,  T. A.  Rupprecht 
201  Lowell  Street 
Wilmington,  Massachusetts 

Raytheon  Company 
Missile  and  Space  Division 
Attn-  Research  Library 
Bedford,  Massachusetts 

Teledyne  Systems,  Incorporated 
Attn  Technical  Library, Antenna  Section 
1625  E.  126th  Street 
Hawthorne,  California 

National  Research  Council 
Attn.  Microwave  Section 
Ottawa  2  Canada 

Sichak  Associates 
Attn:  W,  Sichak 
518  Franklin  Avenue 
Nutley,  New  Jersey 

W,  T .  Patton 
2208  New  Albany  Road 
Cinn  .  Township 
Riverton  Post  Office 
New  Jersey 

Radio  Corporation  of  America 
Missile  and  Service  Radar  Division 
Attn  Manager,  Antenna  Engineering 
Ski  11  Center 
Moore, stown,  New  Jersey 
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Advanced  Development  Laboratories,  Inc. 
Heines  Street 
Nashua,  New  Hampshire 

Cornell  Aeronautical  Lab. 

Attn:  Research  Library 
Buffalo,  New  York 


General  Electric  Company 
Light  Military  Electronics  Department 
French  Road 
Utica,  New  York 

General  Electric  Company 
General  Engineering  Laboratory 
Bldg  371,  Room  478 
Schenectady,  New  York 

Goodyear  Aircraft  Corporation 
Antenna  Department 
Litchfield  Park 
Phoenix,  Arizona 

Laboratory  for  Electronics,  Inc. 
Antenna  Department 
1079  Commonwealth  Avenue 
Boston  15,  Massachusetts 

Lockheed  Aircraft 

Electronic  and  Armaments  System  Office 
Burbank,  California 

Motorola,  Inc. 

Western  Military  Electronics  Division 
P.  0.  Box  1417 
Scottsdale,  Arizona 

Philco  Corporation 
Antenna  Section 
3875  Fabian  Way 
Palo  Alto,  California 


Avco  Corporation 
Electronic  and  Ordnance  Division 
Attn:  Technical  Library 
Cincinnati  41,  Ohio 

Bendix  Corporation 
Research  Labs  Division 
Attn:  G.  M.  Peace 
Southfield  (Detroit),  Michigan 

Douglas  Aircraft  Co.,  Inc. 

Attn:  MSSD  Technical  Library 
(Antenna  Section) 

3000  Ocean  Park  Blvd. 

Santa  Monica,  California 

Emerson  and  Cuming,  Inc. 

Attn:  E.  J.  Luoma 
869  Washington  Street 
Canton,  Massachusetts 


Fairchild  Stratos  Corporation 
Aircraft-Missiles  Division 
Attn:  Technical  Library  (Antenna  Section) 
Hagerstown  10,  Maryland 


Space  Technology  Laboratory 
Attn:  Research  Library 
P.  O.  Box  95100 
Los  Angeles  45,  California 


ANTENNA  LABORATORY 


TECHNICAL  REPORTS  AND  MEMORANDA  ISSUED 


Cont ract  AF33(616)-310 


"Synthesis  of  Aperture  Antennas,"  Technical  Report  No.  1,  C.  T.  A.  Johnk, 
October,  1954.* 

"A  Synthesis  Method  for  Broadband  Antenna  Impedance  Matching  Networks," 
Technical  Report  No.  2,  Nicholas  Yaru,  1  February  1955.  *  AD  61049. 

"The  Assymmetrically  Excited  Spherical  Antenna,"  Technical  Report  No.  3, 

Robert  C.  Hansen,  30  April  1955.* 

"Analysis  of  an  Airborne  Homing  System,"  Technical  Report  No.  4,  Paul  E. 
Mayes,  1  June  1955  (CONFIDENTIAL).  '  ~ 

"Coupling  of  Antenna  Elements  to  a  Circular  Surface  Waveguide,"  Technical 
Report  No.  5,  H.  E.  King  and  R.  H.  DuHamel,  30  June  1955.* 

"input  Impedance  of  A  Spherical  Ferrite  Antenna  with  A  Latitudinal  Current," 
Technical  Report  No.  6,  W.  L.  Weeks,  20  August  1955. 

"Axially  Excited  Surface  Wave  Antennas,"  Technical  Report  No.  7,  D.  E.  Royal, 
10  October  1955.* 

"Homing  Antennas  for  the  F-86F  Aircraft  (450-2500  me),  Technical  Report  No.  8, 
P.  E.  Mayes,  R.  F.  Kyneman,  and  R.  C.  Becker,  20  February  1957,  (CONFIDENTIAL). 

"Ground  Screen  Pattern  Range,"  Technical  Memorandum  No.  1,  Roger  R.  Trapp, 

10  July  1955.*  '  . 

Contract  AF33(616)-3220 


"Effective  Permeability  of  Spheriodal  Shells,"  Technical  Report  No.  9, 

E.  J.  Scott  and  R.  H.  DuHamel,  16  April  1956. 

"An  Analytical  Study  of  Spaced  Loop  ADF  Antenna  Systems,"  Technical  Report 
No.  10,  D.  G.  Berry  and  J.  B.  Kreer,  10  May  1956.  AD  98615. 

"A  Technique  for  Controlling  the  Radiation  from  Dielectric  Rod  Waveguides," 
Technical  Report  No.  11,  J.  W.  Duncan  and  R.  H.  DuHamel,  15  July  1956.* 

"Direction  Characteristics  of  a  U-Shaped  Slot  Antenna,"  Technical  Report 
No.  12,  Richard  C.  Becker,  30  September  1956.** 

"impedance  of  Ferrite  Loop  Antennas,"  Technical  Report  No.  13,  V.  H.  Rumsey 
and  W.  L,  Weeks,  15  October  1956.  AD  119780. 

"Closely  Spaced  Transverse  Slots  in  Rectangular  Waveguide,"  Technical  Report 
No,  14,  Richard  F.  Hyneman,  20  December  1956. 
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"Distributed  Coupling  to  Surface  Wave  Antennas/1  Technical  Report  No.  15 , 

R.  R.  Hodges,  Jr.,  5  January  1957. 

"The  Characteristic  Impedance  of  the  Fin  Antenna  of  Infinite  Length/1 
Technical  Report  No.  16,  Robert  L.  Carrel,  15  January  1957.* 

"On  the  Estimation  of  Ferrite  Loop  Antenna  Impedance/'  Technical  Report  No .  17, 
Walter  L.  Weeks,  10  April  1957„*  AD  143989.  “  "  ”  — 

"A  Note  Concerning  a  Mechanical  Scanning  System  for  a  Flush  Mounted  Line 
Source  Antenna/'  Technical  Report  No.,  18.  Walter  L  Weeks ,  20  April  1957. 

"Broadband  Logarithmically  Periodic  Antenna  Structures , "  Technical  Report  No. 

19 ;  R.  H.  DuHamel  and  D .  E.,  Isbell,  1  May  1957  AD  140734. 

"Frequency  Independent  Antennas/'  Technical  Report  No  20,  V.  H,  Rumsey, 

25  October  1957.  • 

"The  Equiangular  Spiral  Antenna/  Technical  Report  No,  21,  J.  D.  Dyson, 

15  September  1957.  AD  145019.  "  ~ 

"Experimental  Investigation  of  the  Conical  Spiral  Antenna,"  Technical  Report 
No.  22,  R,  L.  Carrel,  25  May  1957-**  AD  144021.. 

"Coupling  between  a  Farallel  Plate  Waveguide  and  a  Surface  Waveguide," 

Technical  Report  No.  23,  E,  J.  Scott,  10  August  1957. 

"Launching  Efficiency  of  Wires  and  Slots  for  a  Dielectric  Rod  Waveguide," 
Technical  Report  No.  24,  J .  W ,  Duncan  and  R,  H.  DuHamel,  August  1957. 

"The  Characteristic  Impedance  of  an  Infinite  Biconical  Antenna  of  Arbitrary 
Cross  Section,1'  Technical  Report  No.  25,  Robert  L,  Carrel,  August  1957. 

"Cavity-Backed  Slot  Antennas,"  Technical  Report  No.  26,  R.  J.  Tector,  30 
October  1957 . 

"Coupled  Waveguide  Excitation  of  Traveling  Wave  Slot  Antennas,"  Technical 
Report  No.  27,  W.  L.  Weeks,  1  December  1957. 

"Phase  Velocities  in  Rectangular  Waveguide  Partially  Filled  with  Dielectric," 
Technical  Report  No.  28 ,  W.  L.  Weeks,  20  December  1957. 

"Measuring  the  Capacitance  per  Unit  Length  of  Biconical  Structures  of  Arbi¬ 
trary  Cross  Section,"  Technical  Report  No.  29,  J.  D.  Dyson,  10  January  1958. 

"Non-Planar  Logarithmically  Periodic  Antenna  Structure",  Technical  Report 
No.  30,  D.  E.  Isbell,  20  February  1958.  AD  156203.  . — 

"Electromagnetic  Fields  in  Rectangular  Slots/1  Technical  Report  No.  31. 

N.  J.  Kuhn  and  P.  E.  Mast,  10  March  1958.  ’  ' 


"The  Efficiency  of  Excitation  of  a  Surface  Wave  on  a  Dielectric  Cylinder," 
Technical  Report  No.  32,  J.  W.  Duncan,  25  May  1958. 

"A  Unidirectional  Equiangular  Spiral  Antenna,"  Technical  Report  No.  33, 

J.  D.  Dyson,  10  July  1958.  AD  201138.  '  '  ^ . . 

"Dielectric  Coated  Spheroidal  Radiators,"  Technical  Report  No!  34,  W.  L. 
Weeks,  12  September  1958.  AD  204547.  — 

"A  Theoretical  Study  of  the  Equiangular  Spiral  Antenna,"  Technical  Report 
No.  35,  P.  E.  Mastj  12  September  1958.  AD  204548. 

Contract  AF33 (616 )-6079 


"Use  of  Coupled  Waveguides  in  a  Traveling  Wave  Scanning  Antenna,"  Technical 
Report  No.  36,  R.  H.  MacPhie,  30  April  1959.  AD  215558. 
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